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Abstract: General dynamic properties like controllability and simulability of 
spin systems, fermionic and bosonic systems are investigated in terms of sym- 
metry. Symmetries may be due to the interaction topology or due to the struc- 
ture and representation of the system and control Hamiltonians. In either case, 
they obviously entail constants of motion. Conversely, the absence of symme- 
try implies irreducibility and provides a convenient necessary condition for full 
controllability much easier to assess than the well-established Lie-algebra rank 
condition. We give a complete lattice of irreducible simple subalgebras of su(2 n ) 
for up to n = 15 qubits. It complements the symmetry condition by allowing for 
easy tests solving homogeneous linear equations to filter irreducible unitary rep- 
resentations of other candidate algebras of classical type as well as of exceptional 
types. — The lattice of irreducible simple subalgebras given also determines mu- 
tual simulability of dynamic systems of spin or fermionic or bosonic nature. We 
illustrate how controlled quadratic fermionic (and bosonic) systems can be sim- 
ulated by spin systems and in certain cases also vice versa. 
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1. Introduction 

Experimental control over quantum dynamics of manageable systems is par- 
amount to exploiting the great potential of quantum systems. Both in simu- 
lation and computation the complexity of a problem may reduce upon going 
from a classical to a quantum setting [H212] ■ On the computational end, where 
quantum algorithms efficiently solving hidden subgroup problems U have es- 
tablished themselves, the demands for accuracy ('error-correction threshold') 
may seem daunting at the moment. In contrast, the quantum simulation end 
is by far less sensitive. Thus simulating quantum systems [5] — in particular at 
phase-transitions [HIE] — has recently shifted into focus [51I WTUl[TTl[r2] . In view 
of experimental progress in cold atoms in optical lattice potentials (T5irPf] as 
well as in trapped ions |15U16j , Kraus et al. have explored whether target quan- 
tum systems can be universally simulated on translationally invariant lattices of 
bosonic, fcrmionic, and spin systems [17] . In some respect, their work can also 
be seen as a follow-up on a study by Schirmer et al. [18] (see also recent work by 
Wang et al. [TO]) specifically addressing controllability of systems with degen- 
erate transition frequencies. Many experimental tasks are engineering problems 
that profit from quantum systems theory as a framework and optimal control 
algorithms for solving the actual problem. 

As compared to an abstract point of view |20j . the flavour of quantum sys- 
tems theory pursued here is meant to be very pragmatic: it takes the causal 
formulation of dynamic systems |21) and does not care about specifics of the 
quantum measurement problem beyond the basic notions [22] and some recent 
developments [23] . Yet it is for these reasons that quantum systems and con- 
trol has quite generally been recognised as a key generic tool [24 ,25,26] needed 
for advances in experimentally exploiting quantum systems for simulation or 
computation and even more so in future quantum technology. It paves the way 
for constructively optimising strategies for experimental implcmentions in real- 
istic settings. Moreover, since such realistic quantum systems are mostly beyond 
analytical tractability, numerical methods are often indispensable. To this end, 
gradient flows can be implemented on the control amplitudes thus iterating an 
initial guess into an optimised pulse scheme 27,28,29 . This approach has proven 
useful in spin systems [20] as well as in solid-state systems [21] ■ Moreover, it has 
recently been generalised from closed systems to open ones [32) , which are known 
to be a challenge to control [53], where the Markovian setting can also be used 
as embedding of explicitly non-Markovian subsystems [34] . 

However, in closed systems, the numerical tools usually require the system is 
universal or fully operator controllable [3"5ll3"6] . For a plethora of systems with 
symmetry constraints we have recently determined explicit dynamic system alge- 
bras [23 (as subalgebras of su(N)), and conversely, we have derived design rules 
for the experimenter as guidelines ensuring universality of quantum architec- 
ture. While extending earlier work on branching diagrams of simple subalgebras 
of su(N) S38ll39j . here we focus on complete necessary and sufficient conditions 
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for full controllability (mostly) confining ourselves to arguments easy to check 
by inspection or to decide by computationally cheap algorithms such as solving 
a system of homogeneous linear equations. 

In view of applications, we illustrate our findings by a comprehensive set of 
worked examples on spin chains. Actually Ising-ZZ coupled n-spin-i chains with 
mostly collective controls or Hciscnbcrg-XX chains with one single local control 
suffice to get exponential growth of dynamic degrees of freedom (in the sense 
their respective dynamic system algebras are sp(2 n_1 ) or so(2™)). Our work thus 
adds to the recent spin-chain literature (see, e.g., [i01HTlH2lH51l441l45U46HTnH51 
ITT)] and compare [49,50,51]) and — on a more general scale — it is anticipated to 
have significant impact on quantum simulation as well as distributed quantum 
computing (see, e.g., [5"2"Il5"3T5"4"] ) . 

2. Overview and Main Results 

More precisely, the first main part develops, starting from the basic notions of 
controllability (Sec. [3]) in terms of coupling graphs (Sec.HJ) and their symmetries 
(Sec. [5]), a single necessary and sufficient symmetry condition for full control- 
lability (Sec. [7]). To this end and in view of practical applications, Sec. [6] gives 
branching diagrams of all irreducible simple subalgebras of the unitary algebras 
su(N) with N < 2 15 . Concomitantly, we provide a set of efficient computational 
algorithms for assessing controllability by merely solving systems of homoge- 
neous linear equations. 

The second part focusses on simulability (Sec. [5]) in terms of dynamic sys- 
tem algebras. A plethora of worked examples is discussed in Sec. including 
four full series of qubit chains coupled by pair interactions such that their dy- 
namic system algebras for the first three cases are so(2n + I), so(2n + 2), and 
sp(2™ _1 ), respectively. Most remarkably, for n > 4, the fourth series results in 
dynamic system algebras so(2 n ) if (n mod 4) G {0,1} and sp(2 n ~ 1 ) else. The 
findings also interrelate spin systems, fermionic systems (Sec. ITU]) and bosonic 
systems (Sec. [TTJ). The algebraic conditions for simulability given are sufficient 
to ensure the existence of solutions to the actual task of quantum simulation of 
closed systems formulated as an observed optimal control problem in the outlook 
(Sec.UJ). 

3. Controllability 

Consider the controlled Schrddinger equation lifted to unitary maps (quantum 
gates) 

m 

U(t) = -t(H d + J2'u0)Hj)U(t) . (1) 
j'=i 

Here the system Hamiltonian Hd denotes a non-switchable drift term and the 
control Hamiltonians Hj can be steered by (piece-wise constant) control am- 
plitudes Uj(t) € K taken to be unbounded henceforth. The equation of motion 
governs the evolution of a unitary map of an entire basis set of vectors represent- 
ing pure states. Using the short-hand notations H := Hd + J2]Li u j(t)Hj an d 
adir(vec/l) := [H,A], the Liouville equation p{t) = —i[H,p(t)] can be rewritten 



vecp(i) = — i adff vecp(t). 



(2) 
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Algorithm 1: Determine system algebra via Lie closure 

Input: Hamiltonians / := {iH^; iHi, . . . , iH m } C su(TV) 

1. B := maximal linearly independent subset of I 

2. num := #B 

3. If num = TV then O := B else O := { } 

4. If num = TV or #B = then terminate 

5. C := [O, B] U [B, B], where 
[Si,S 2 ] = {[s 1 ,s 2 ]\s 1 6Si,s 2 €5 2 } 

6. O := OUB 

7. -B := max. linear independent extension of O 
with elements from C 

8. num := num + #B; Go to 4 

Output: basis O of the generated Lie algebra and 
its dimension num 

The complexity is roug hly 0{N 6 ■ TV 2 ), as about TV 2 times 
a rank-revcaling QR decomposition has to be performed in 
Liouville space (with dimension TV 2 ). For n qubits, TV := 2™. 



Both equations of motion take the form of a standard bilinear control system 
(£) known in classical systems and control theory |55) 

m 

X(t)=(A + J2u j (t)B j )X(t) (3) 
i=i 

with 'state' X(t) e C^, drift A e qI(N,C), controls B 3 £ g[(iV,C), and control 
amplitudes Uj £ K, where gl(N, C) denotes the set of complex TV x TV matrices. 
Since all the control systems considered henceforth are bilinear, we often drop the 
specification bilinear for short. Now lifting the (bilinear) control system (£) to 
group manifolds [551157] by X(t) E GL(tV, C), i.e. the set of non-singular complex 
TV x TV matrices, under the action of a compact connected Lie group K with Lie 
algebra % while keeping A, Bj € g[(TV, C), the condition for full controllability 
turns into the Lie algebra rank condition |58U591l57] 

(A^|i = l,2,...,m) L ie = «, (4) 

where (-)Lic denotes (the linear span over) the Lie closure obtained by repeatedly 
taking mutual commutator brackets. Algorithm 1 gives an explicit method to 
compute the Lie closure, see also [50] . 

Transferring the classical result [555] to the quantum domain [5T1I^1[5d] . the 
bilinear system of Eqn. (JTJ) is fully (operator) controllable if and only if the drift 
and controls are a generating set of the special unitary algebra su(N): 

(iH d , iHj | j = 1,2,..., m) Li c = « = su(iV). (5) 

In fully controllable systems, to every initial state po the reachable set is the 
entire unitary orbit Ou(po) '■— {U p^W \ U S SU(tV)}. With density operators 
being Hermitian this means any final state p(t) can be reached from any initial 
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state po as long as both of them share the same spectrum of eigenvalues. Thus 
reachable sets and isospectral sets coincide. 

In contrast, in systems with restricted controllability the Hamiltonians gen- 
erate but a proper subalgebra of the full unitary algebra 

(iH d , iHj | j = 1,2,..., m) Uc = I C su(N). (6) 

Then the dynamic group K := expt is but a proper subgroup K C SU(iV) of 
the full unitary group. Therefore the corresponding reachable sets take the form 
of subgroup orbits of initial states 

teacl)Oo) = OkM ■= {Kp Q K^ \KeKC SU(JV)}. (7) 

4. Natural Tensor-Product Structure and Coupling Graphs in Qubit 
Systems with Pair Interactions 

We start out with the case of qubit systems coupled by pair interactions. Yet 
quantum simulation of effective many-body interactions in multi-level systems 
requires more refined notions, see Appendix lAl and 151 Thus we choose a line-of- 
thought allowing for the extensions needed later in a natural way while trying 
to keep the overhead minimal here. Finally it should be stressed that the results 
in Sees. [5H7] are valid in full generality of Appendix lAl and [B] 

To fix the basic terminology, observe that the abstract direct sum of Lie alge- 
bras has a matrix representation as the Kronecker sum, e.g., su(di) © su(<i2) •= 
su(e?i) <g) td 2 + Id! <8> su(e?2) and that it generates a group isomorphic to the 
Kronecker product (i.e. tensor product) G = SU(di) <g> SU^). The abstract 
direct sum of two algebras f)i and f)2 (each given in an irreducible representa- 
tion) has itself an irreducible representation as a single Kronecker sum f)i ffi ()2 
(Thm. 11. 6. II of Ref. [63]). Such an irreducible direct sum representation always 
exists for every semi-simple Lie algebra which is not simple. 

Control systems consisting of n qubits are usually embedded in su(N) with 
N := 2". Their natural intrinsic tensor-product structure takes the form of the 
n-fold Kronecker sum su(2)©su(2)© ••• ©su(2). An A^ 2 -l dimensional skew- 
Hermitian tensor basis with respect to this tensor-product structure can be given 
in terms of the Pauli matrices 

1 := h = (o i) ' x := a * ^ (? o) ' Y := °y = (S 7) ' and z := a * = (0 -1) 

(8) 

by defining the elements — 5H1H2 • • • H n , where HiH 2 • • • H„ := Hi®H2®- • -(8H n 
and Hj £ {I, X, Y, Z}. The element Hi = H2 = ■ ■ ■ = H„ — I is not traceless and 
hence cannot occur in su(2"). In terms of this tensor basis, we write Hamiltonians 
as linear combinations (c^ £ R) 

m 

H = J2c k U k (9) 

fc=i 

of elements U k = -\i^ih,x®'rih,%®- • ■®'Hk,n) with^fe^ £ {I,X,Y, Z}. Consider- 
ing local controls and pairwisc coupling interactions the orders of the constituents 
are confined, i.e. 



ord(H fe ) := #{£ : U k ,e + h} £ {1,2}. 
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Fig. 1. General coupling topology represented by a connected graph. The vertices denote the 
spin-i qubits, while the edges represent pairwise couplings (e.g. of Heisenberg or Ising type). 
Qubits of the same colour and letter are taken to be affected by joint local unitary operations 
(or none if the color is white), while qubits of different kind can be controlled independently. 
For a system to show an outer symmetry brought about by permutations within subsets of 
qubits of the same type, both the graph as well as the system plus all control Hamiltonians 
have to remain invariant. 



Usually, the control Hamiltonians Hj are local, i.e. all terms in Eqn. © (for 
H = Hj) arc of order one, while the corresponding terms in Eqn. (|9|) for the 
drift Hamiltonian Hq (= H) arc of order two comprising the non-switchable 
pairwise coupling terms. 

Now, in a coupling graph the vertices representing the local subsystems are 
connected by edges, where each edge stands for a pairwise coupling term occuring 
in the drift Hamiltonian Hd- An example of a connected coupling graph is shown 
in Figure [TJ — Connected coupling graphs are essential for full controllability as 
elucidated by the following theorem. 

Theorem 1. Consider a bilinear control system with pair interactions on su(2 ra ), 
where all the local subsystems su(2) are independently fully controllable so the 
dynamic algebra I D su(2) ©su(2) © ■ • ■ ©su(2). Then the system is fully control- 
lable, i.e. t = su(2 n ), if and only if its coupling graph is connected. In particular, 
t = su(2 n ) is simple. 

Proof. A proof is given in Ref. [64] (see Thm. 2, Remark 5.1, and Thm. 4), see 
also Ref. [52]. ■ 



5. Symmetry-Constrained Controllability 

A Hamiltonian quantum system is said to have a symmetry expressed by the 
skew-Hcrmitian symmetry operator s € Su(N), if 

[s,H v ]=0 for all v £ {d; 1,2,..., m}. (10) 

More precisely, we use the term outer symmetry if s generates a SWAP operation 
permuting a subset of qubits of the same type (cp. Fig. [TJ) such that the coupling 
graph and all Hamiltonians {H v } are left invariant. Now subsets of qubits are 
termed indistinguishable if and only if they can be interchanged by an outer 
symmetry, i.e. a SWAP operation that is a symmetry of the system; otherwise 
they are distinguishable. In contrast, an inner symmetry relates to elements s not 
generating a SWAP operation in the symmetric group of all qubit permutations. 
In either case, a symmetry operator is an element of the centraliser 



{H v }' := Z BU(N) ({H„}) = {s e su(N) \ [s, H v ] = \/v G {d; 1,2, . . . , to}}, 

(11) 
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Algorithm 2: Determine centraliser (resp. commutant) 
to system algebra t 

Input: Hamiltonians / := {iH^ \ iH\ , . . . , iH m } C M 

1. For each Hal solve the homogeneous linear eqn. 
S H := {s eM\(l®H-H*® 1) vccO) = 0} 

2. _R := f] HeI S H . 
Output: 

R _ f centraliser I' if M = su(N) 

~ 1 commutant of t if M = gl(N, C) 



The complexity is roughly 0(N 6 ), as in Liouville space 
N 2 equations have to be solved by LU decomposition. 
For n qubits, N := 2". 



recalling that the centraliser of a given subset m C su(7V) with respect to a Lie 
algebra su(JV) consists of all elements in su(N) that commute with all elements 
in m. Jacobi's identity [[a, b], s] + [[b, s], a] + \[s, a], b\ —Q gives two useful facts: 
(1) an element s that commutes with the Hamiltonians {iH v } also commutes 
with their Lie closure t. For the dynamic Lie algebra t we have 

«' := Z m(N) {l) = {s G Su(iV) | [a, k] =0 Vfc £ t} (12) 

and hence {iH„}' = t' . Thus in practice it is (most) convenient to just evalu- 
ate the centraliser for a (minimal) generating set {iH^} of t since the overall 
symmetry properties can be read from the local symmetries of the constituent 
Hamiltonians. Fact (2) means the centraliser 6' forms itself an invariant Lie sub- 
algebra (or ideal) to su(iV) collecting all symmetries. In summary, we obtain the 
following straightforward, yet important result: 

Theorem 2. Lack of symmetry in the sense of a trivial centraliser is a necessary 
condition for full controllability. 

Proof. Any non-trivial element in the centraliser would generate a one-parameter 
subgroup in K' c SU(iV) that is not in K = cxpt. ■ 

Throughout this paper, we consider finite-dimensional complex matrix repre- 
sentations of Lie algebras, a representation being a map from a given Lie algebra 
to the set complex square matrices of appropriate (and finite) dimension. The 
matrix entries are given by complex polynomial (or equivalently holomorphic) 
functions. In the following, we will usually not consider the trivial representa- 
tion, which maps any element to 1 £ C. One particular important example for 
a representation of a Lie algebra is the standard representation, which is the 
lowest-dimensional (non-trivial) representation (with some exceptions, see the 
Appendix [C| and which is typically used to define the corresponding Lie algebra 
in its matrix form. In analogy to the centraliser, one can define the commutant 
relative to a representation <j> of dimension dim(0) 



comm^ (m) := {g € g [(dim ((/)), C) | \g,(j)(m)] = Vm £ m} (13) 
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for a subset m C g of a Lie algebra g. Now it is natural to ask how the notions 
of centraliser and commutant relate to irreducible representations. 

Lemma 3. Let <S> denote the standard representation of su(N). If 6 C su(N), 
then the following statements are equivalent: 

(1) The centraliser 6' = -Z su (JV)(t) oft in su(N) is trivial, i.e. zero. 

(2) The restriction of <P from su(N) to 6 is irreducible. 

(3) The commutant comm$(f) of t w.r.t. <P is trivial, i.e. = {c ■ In |c £ C}. 

Proof. As su(N) is compact, it follows that <P and its restriction to 6 are com- 
pletely reducible in the sense of being a direct sum of irreducible representa- 
tions (see Cor. 2.17 of [65]). The representation <£> is even irreducible and faith- 
ful, i.e. injective. Hereafter, we will consider the complcxification 6c of 6 and 
su(N)c = si(N,C) as complcxification of su(N). The representation <P has a 
unique extension <t>c to sl(N,C), which is also irreducible and faithful. In addi- 
tion, <l>c and its restriction to 6c are completely reducible. These facts can be 
deduced from Thm. 1, pp. 111-112 of [55] and Prop. 7.5 of |6"T] . 

Now it follows that (1) is equivalent to Z s i(N)(tc) = {0}- As @c is faithful, 
this holds if and only if comm# c (6c) is trivial. Relying on the fact that <?c is 
completely reducible, comm<£ c (6c) is trivial if and only if the restriction of <?c 
from sl(N,C) to 6 C is irreducible. Using Thm. 1, pp. 111-112 of [5J], this is 
equivalent to (2). As $ is completely reducible, (2) and (3) arc equivalent. ■ 

As a second consequence of a trivial centraliser the corresponding subalgcbra 
6 of su(N) has to be simple or semi-simple: 

Lemma 4. Let 6 C su(N) be a subalgebra to the Lie algebra su(N). If its cen- 
traliser 6' in su(N) is trivial, then 6 is simple or semi-simple. 

Proof. By compactness, 6 = 3e © s decomposes into its centre it and a semi- 
simple part 5 (see, e.g., Cor. IV. 4. 25 of Ref. [67]). As the centre it = 6' PI 6 is 
trivial, 6 can only be semi-simple or simple. M 

Note that the centraliser is 'exponentially' easier to come by than the Lie 
closure in the sense of comparing the asymptotic complexity 0(N 6 ■ N 2 ) (with 
N := 2™ for n qubits) of Algorithm 1 for the Lie closure with the asymp- 
totic complexity 0(N 6 ) of Algorithm 2 for the centraliser tabulated above. — 
Therefore one would like to fill the gap between lack of symmetry as a neces- 
sary condition and sufficient conditions for full controllability in systems with 
a connected coupling topology. For pure-state controllability, this was analysed 
in [64) . for operator controllability the issue has been raised in [25], inter alia 
following the lines of |681l69j , however, without a full answer. Further results in 
the case of pure-state controllability can be found in [33] . 

We have proven that the lack of symmetry is necessary for a control system 
to be fully controllable. Yet in turn, a control system without symmetry need 
not be fully controllable, as the following elementary (and pathological) example 
shows: 

Example 5. Assume we have a bilinear control system on two qubits, where the 
dynamic Lie algebra 6 = (iXI, iYl, iZI, ilX, ilY, iIZ)uc = su(2)0su(2) is not 
simple. Although it has no symmetry and its centraliser 6' is in fact trivial, the 
system is not fully controllable: all pair terms like iZZ cannot be generated, since 
its pathological 'coupling graph' 
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is clearly not connected. 

Nevertheless, the somewhat trivial example is illuminating. While in the con- 
text of C*-algebras, von Neumann's doublc-commutant theorem recovers the 
original algebra from the commutant of its commutant |70U71j , a similar theorem 
does not extend to Lie algebras [72]. Rather, if the dynamic algebra t C su(N) 
has a trivial centraliser t' , then the double centraliser t" , i.e. the centraliser of the 
centraliser in su(N), of all compact semi-simple and simple irreducible proper 
and improper subalgcbras I of su(N) is given by su(iV) in line with Lemma 2J 
However, if one considers the associative matrix algebra (with identity) gener- 
ated by the basis elements (including the identity matrix) of a Lie algebra via 
its standard representation, then von Neumann's double commutant theorem 
still holds, see Thm. (3.5.D) of Ref. [75]. — In the next step, we will thus add a 
criterion to single out the simple subalgebras. 

Motivated by Example [5] one might conjecture that the dynamic algebra t 
is simple if £ acts irreducibly and the coupling graph of the control system is 
connected. This is true for control systems in qubits with pairwise coupling 
interactions: 

Theorem 6. Consider a bilinear control system with pair interactions onsu(2 n ). 
Assume that the tensor-product structure is given by su(2) ©su(2) © • • • ffisu(2) 
and that the centraliser t' of the dynamic algebra t is trivial. The dynamic algebra 
t is simple if and only if the coupling graph of the control system is connected. 

Proof. See Corollary 2) in Appendix [B] ■ 

The general case beyond pair interactions (and qubit systems) is discussed in 
Appcndix[B] In the case of pair interactions, we say a control system is connected 
if its coupling graph is connected. This definition of a connected control system is 
a particular case of the general definition (see Appendix [Bj applicable to control 
systems which do not have a natural coupling graph. 

6. Irreducible Simple Subalgebras of su(N) 

Starting from the knowledge that for a fully controllable system the dynamic 
algebra t has to be simple and given in an irreducible representation (see, e.g., 
Appendix IB"]) . it is natural to ask for a classification of all these cases. Following 
the work of Killing, Elie Cartan [73] classified all simple (complex) Lie alge- 
bras (see, e.g., [751175] ). The corresponding compact real forms ( [751177] ) are the 
compact simple Lie algebras of classical type (assuming £ £ N \ {0} henceforth): 

ae : su(£ + 1), 
b e : so{2t + l), 

c e : sp{£) :=sp(2^,C) Hu(2^,C), 
Of : So(2£), 

and of exceptional type t^, tr, t8: ?4- 02- Note also that for ae (£ > 1), be (£ > 2), 
ce (£ > 3) and Ve > 4) the following isomorphisms (see, e.g., Thm. X.3.12 in 
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[77] ) su(2) = so(3) = sp(l), so(5) = sp(2), and su(4) = so(6) arc no longer of 
concern. The same holds for the abelian case 50 (2) as well as for the semi-simple 
oneso(4) S su(2) ©su(2). 



Pro memoria. The classical, compact simple Lie algebras and some forms of their 
standard matrix representations may be given as follows: 



algebra 


definition and block forms 


Lie dimension 


su(iV) := {a e 


C NxN | t tra = } 


iV 2 -l 


so(N) := {a = 


Uatf e C NxN | 5* = -a, tra = 0} 
with U e SE/(./V) 


|JV(JV — 1) 


N = 21: 


, (A B\ 
a = [c -A 1 ) 






with A, B, C e C £x£ , = -5, C* = 


-c 


N = 21+1: 


a' = C -A 1 v\ 

\-it* -u* oy 

with A, £>, C as above and u, v £ C 




NB: the representations a and a' above need not be equal. 





sp(7V/2) := {a = UaU^ e C NxN | Ja = -a* J} ±JV(JV + 1) 

with J:=(^ -**),U€SU(N) 



with 4,B,Ce C^, B = B\C = C t 



By completeness 1 of Cartan's classification above we may summarise as follows: 

Corollary 7 (Candidate List). Consider a bilinear control system, where the 
drift and control Hamiltonians {iH„} generate the dynamic system Lie algebra 
6 C su(N) in an irreducible representation (f trivial) with the additional promise 
that t is simple (e.g. due to a connected control system). Then, being a simple 
subalgebra ofsu(N), the system algebra t has to be one of the candidate compact 
simple Lie algebras: su(£ + 1), so(2£ + 1), sp(£), so(2£), t 6 , tf, eg, f4, or q 2 - • 

For illustration of the Lie algebras of exceptional type, consider the dimensions of 
their standard representations (see, e.g., p. 218 of Ref. [75], Ref. [73], Ref. [50] ) 
e 6 C su(27), e 7 C sp(28), e 8 C so(248), f 4 C so(26), and q 2 C so(7). As a 
final remark on exceptional Lie algebras suffice it to add that — with the single 



1 NB: The list of algebras is indeed complete - note that in particular spin and pin groups 
are also generated by the algebras bo(N) and o(N), respectively 1781 . 
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Table 1. The Irreducible Simple Subalgebras of su(JV) for N < 16 



eu(2) 



su(2 



su(2)«-sp(2 



su(2)«-so(5 



su(2)«-sp(3 
so (6 
su(3 



su(2)«- 02 <-so(7 



su(2)«-sp(4 

su(3)*, 

\-so 8 
so (7) J 



su(2)<-so(9 



-su(3) 



-su(4) 



-su(5) 



■su(6) 



-su(7) 



■au(8) 



-su(9) 



su(2)«— sp(5)«i 
so(5)<-so(10)<- 
su(3)«-- 
su(4)<- 
su(5)<- 



su(10) 



su(2)<-so(ll) 


«-*u(ll) 


su(2)—sp(6)< 


Jjtt(12) 




S0(12) 


6U(2)«-*>(13) 


<-su(13) 



sp 3 J 



so(5) 
sp(3) 
g2«- 



■so(14)' 



■5u(14) 



kso(15V 
so(6)J 

su(3) + 
su(5)<- 
5u(3) «— su(6) «- 



■su(15) 



|-sp(8)«-, 



sp(2)< 

so(9)«-so(16)«|" su ( 16 ) 
so(10)« 



NB: so(3) su(2) ^ sp(l), so(5) sp(2), and so(6) S su(4). 



exception of Q2 — they all fail to generate groups acting transitively on the sphere 
or on R-^ \ {0}. This has been shown in [5T] building upon results in [52] to fill 
gaps in earlier work [551151] . 

Having listed all the candidates for proper simple subalgebras of su(N), wc 
now focus on the set of possible irreducible representations. To this end, in this 
chapter we describe the main results, while all the details shall be explained in 
the Appendix [C] The irreducible representations of simple (complex) Lie alge- 
bras were already determined by Elie Cartan |85j . This classification is equiva- 
lent for the compact simple Lie algebras (or the compact connected simple Lie 
groups), see, e.g., [75]. The irreducible simple subalgebras of su(N) are found 
by enumerating for all simple Lie algebras all their irreducible representations of 
dimension N. The dimensions of the irreducible representations can be efficiently 
computed using computer algebra systems such as LiE [55] and MAGMA [57] via 
Weyl's dimension formula. Following the work of Dynkin [55] (see App. IC.3l and 
Chap. 6, Sec. 3.2 of Ref. [HO]); one can determine the inclusion relations between 
irreducible simple subalgebras of su(N) . We obtained all the irreducible simple 
subalgebras of su(N) for N < 2 15 = 32768. This significantly extends previous 
work [551155] for N < 9. The results for N < 16 are given in Tab. [1] those for 
N = 2™ and 1 < n < 15 are relegated to Tab. [2] A complete list with all the 
results for N < 2 15 is attached as Supplementary Material [91] , 

Remark 8. With regard to Tables [T] and [5] the occurrence of su(2) as an ir- 
reducible simple subalgcbra to any su(iV) with N > 2 is natural from the 
point of view of spin physics. Wc identify su(N) = su(2j + 1), where the (non- 
vanishing) half-integer and integer spin-quantum numbers may take the values 
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j e {i, 1,|,2,...}. Now to any such j there is an irreducible spin-j represen- 
tation of the three Pauli matrices generating su(2). For instance, in su(4) there 
is an irreducible spin-| representation of su(2) as a proper irreducible subalge- 
bra su(2) C su(4). — In contrast, the Gell-Mann basis to su(2j + 1) comprises 
a reducible representation of su(2) as a subalgebra. Clearly, the two types of 
representations are mequivalcnt. 

In the set of irreducible simple subalgcbras of su(N), the subalgcbras sp(N/2) 
with N even and so (AT) play a particularly important role. For N > 5, we discuss 
the irreducible simple subalgcbras of su( N) for N even and odd. If N > 5 is even, 
then su(N) has both sp(N/2) and so(N) as irreducible simple subalgebras. In 
addition, su(2) C sp(AT/2) occurs as irreducible simple subalgebra. We consider 
two types of trivial cass. First, if N > 5 is even and if sp(N/2), so (AT), and 
su(2) C $p(N/2) are the only proper irreducible simple subalgebras, then we 
say the case is trivial. A trivial example is given by su(12) in Tab. [TJ If A^ > 5 
is odd, then so (AT) is an irreducible simple subalgebra of su(N) but sp(AT/2) is 
not (as A^/2 is not an integer). Moreover, su(2) C so(N) occurs as irreducible 
simple subalgebra. Second, if A" > 5 is odd and if so (AT) as well as su(2) C 
so (AT) are the only proper irreducible simple subalgebras, then we say the case 
is trivial. Examples of such trivial cases are given by su(5), su(9), su(ll), and 
su(13) in Tab. [TJ The irreducible subalgebras sp(N/2) and so(A^) correspond 
to the symmetric spaces SU(AT)/Sp(Ar/2) and SU(AT)/SO(Ar). These are two of 
three possible symmetric spaces [77] of SU(A r ), where the third type docs not 
correspond to a semi-simple subalgebra of su(N). 

We call a representation cf> of a subalgebra 6 symplectic [orthogonal] if the 
subalgebra 6 given in the representation <j> is conjugate to a subalgebra of sp(N/2) 
[so(N)}. If the representation is neither symplectic nor orthogonal, we term it 
unitary. In abuse of notation, we call also the subalgebra i (w.r.t. some fixed but 
unspecified representation (f>) symplectic, orthogonal, or unitary, if the respective 
representation (j> is symplectic, orthogonal, or unitary. 2 We emphasise that the 
classification of a subalgebra depends on the representations considered, sec also 
Chap. IX, App. II.2, Prop. 3 of Ref. [75]. 

The property of a representation to be symplectic [resp. orthogonal] corre- 
sponds to the existence of an invariant (non-degenerate) skew-symmetric [resp. 
symmetric] bilinear form on the space of the representation. For irreducible repre- 
sentations in the compact case [e.g. for subgroups of SU(A f )], this correspondence 
is an equivalence and a proof can be found in Sec. 3.11, Thm. H of Ref. [55] . 
As an invariant (non-degenerate) bilinear form can either be skew-symmetric or 
symmetric, it follows that the same holds for the classification of (irreducible) 
symplectic, orthogonal, or unitary representations (Chap. IX, Sec. 7.2, Prop. 1 
of Ref. [75]): 

Lemma 9. An irreducible representation (/>(?) can either be symplectic, or or- 
thogonal, or unitary. 



2 This notation is motivated by the classification of representations and subalgebras as sym- 
plectic and orthogonal in Ref. 1921 and in Chap. VIII, Sec. 7.5, Def. 2 of Ref. 1761 . Classifying 
representations and subalgebras as unitary appears to be non-standard notation. Unfortu- 
nately, the respective representations are also said to be of quaternionic, real, or complex 
type. 



Symmetry Principles in Quantum Systems Theory 

Table 2. The Irreducible Simple Subalgebras of su(2 n ) for 1 < n < 15 



13 



eu(2) 



su(2)<-sp(2)<-su(4) 


su(2)<-sp(4)<- 




su(3)^ 


-su(8) 


[to 8 )J 




so (7 J 





J-sp(8)«n 



*U(2)«-| 
sp(2) 
so(9)<-so(16) 
so(10) 



■su(16) 



5u(2)<-| 

r-*P(16)«i 
5o(ll)«-so(12)J Uu(32) 
so(32)J 




<su(2)«— sp(64)«-i 
so(9)<-, ksu(128) 
^so(16)«-so(128)* J 

80(15)^ 

«u(2)<-i 

k«n 128 )« 
ep(2)«-l 

ao(17)- S o(256)j su(256) 
5u(4)<- 
so(18)< 



«u(2) 
so(19)«-so(20) 



1-sp(256) 



bu(3) 

sd(7)<-|-bo(512)- 
sp(3)< 



■su(512) 



su(2)- 

5o(21)<--sp(512)» 
sp(2)J 

su(5)<-so(1024) 
5o(22)«- 



i3u(1024) 



su(2)<— sp(1024) 
so(23)<-so(24)«-so(2048) 



j-su(2048) 



Hsp(2048) 



su(2)<-, 

sp(2). 

su(3)* 

8o(7> 
so(17) 
6o(25) 

sp (4)J-so(4096) 



60 (6) 

so(8)<- 

02" 



60(26) « 



£su(4096) 



6u(2)«-| 

^sp(4096)i 
so(27)<-so(28)J h u ( 819 2) 
so(8192)*J 



su(2). 
so(29) T sp(8192)< 



sp(2). 



so(16384)- 
su(4)«— 
5o(30)<- 



•su(16384) 



eu(2) 
su(6) 



su(3) 
sp(3)<- 
8o(8)<- 
;so(31)<-so(32) 



J-sp(16384) 



■so(32768)- 



au(32768) 
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7. Prom Necessary to Sufficient Conditions for Controllability 

While the ramification of mathematically admissible irreducible simple candidate 
subalgebras may seem daunting, in the following we will eliminate candidates 
by simple means. More precisely, we arrive at the following. 

Corollary 10 (Task List). One way of showing full controllability amounts to 
excluding other candidates of irreducible simple subalgebras, which can be 

(1) symplectic, i.e. conjugate to a subalgebra of sp(N/2), 

(2) orthogonal, i.e. conjugate to a subalgebra ofso(N), 

(3) or unitary in the remaining cases. 

In particular, one has to exclude cases like the exceptional ones cq, z-j, es f4, Q2- 
The unitary, irreducible simple subalgebras can occur in the cases su(£ + 1) C 
su(JV) (t>2), so(U + 2), andc 6 . ■ 

In what follows, the plan is to make use of the fact that in Tabs. [T] and [H most 
of the irreducible subalgebras are symplectic or orthogonal. The symplectic and 
orthogonal ones (including their nested subalgebras!) will be excluded by merely 
solving simultaneous systems of linear homogeneous equations, which will also 
exclude the exceptional algebras e?, es, fi, and 02, just leaving Cq. It appears 
that for systems of dimension 2", irreducible representations of cannot be 
an irreducible simple subalgebra of su(2 n ) without being a subalgebra to an 
intermediate orthogonal or symplectic algebra. 

In principle, the task of identifying the dynamic Lie algebra can also be solved 
by algorithms [Mll§5] available in the computer algebra system MAGMA [57]. Yet, 
here we focus on exploiting algorithms that are more efficient, since they boil 
down to solving systems of homogeneous linear equations, which currently can — 
in general — be carried to matrix sizes of about 6 • 10 4 x 6 • 10 4 (and in extreme 
cases to 10 5 x 10 5 )[95]. So the algorithms presented here aim at the more specific 
task of distinguishing su(N) from its proper irreducible simple subalgebras, a 
task our algorithms arc more efficient in. 



7.1. Symplectic and Orthogonal Subalgebras. In order to decide on conjugation 
to irreducible subalgebras which are symplectic and orthogonal, we need more 
detail. To this end 3 , recalling the following Lemma will prove useful to apply the 
lines of |99j in streamlined form leading to an explicit algorithm. 

Lemma 11. (1) Every unitary symmetric matrix S = S l £ U(N,C) is unitarily 
t-congruent to the identity, i.e. S — T IT with T unitary. 

(2) Every unitary skew- symmetric matrix S = —S G U(N,C) with N even is 
unitarily t-congruent to J, i.e. S = T l JT with T unitary and 

MC'H (I4) 

Proof. (1) Follows by singular-value decomposition and goes back to Hua QlOOj . 
Thm. 5). (2) Follows likewise from the same source (ibid., Thm. 7). ■ 



3 Preliminary results were given in the conference papers 97 , 98 . 
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Algorithm 3: Check conjugation to subalgebras of so(N) or sp(JV/2) 



Input: Hamiltonians I := {iH^; iH%, . . . , iH m } C su(N) 

1. For each Hamiltonian Hal determine all non-singular 
solutions to the homogeneous linear equation 

S H ■= {S G SL(N)\SH + H l S = 0} 

= {S G 5L(Af)|(H* <gi 1 + 1 ® vec(S) = 0} 

2. 5 := fine/ 5 » 

Ouiput; (a) 3S G 5 s.t. 5S = +1 <^ t C so(7V) 

(b) 35 G 5 s.t. SS = -1 <=> f C sp(N/2) 

(c) ^S G 5 => t £ so(AT) and { g sp(N/2) 
The cases (a) and (b) are mutually exclusive 
if the centraliser of / is trivial. 



The complexity is roughly 0(7V 6 ), as in Liouville space 

N 2 equations have to be solved by LU decomposition (N := 2 n ). 



Lemma 12. Suppose t C Su(N) is simple and J is defined as in Eqn. (|14[) . Then 
the element iH G t 

(1) is unitarily conjugate to iH e so(N), where H = —H, if and only if there 
exists a symmetric unitary S (so SS = +1jv ) satisfying SH + H S = 0; 

(2) is unitarily conjugate to iH € Sp(N/2) (with N even), where JH = — H l J , 
if and only if there is a skew- symmetric unitary S (so SS = — In) satisfying 
SH + H l S = 0. 

Proof. First observe that whenever there is a unitary T such that THT^ =: H 
with LH = —H l L, this is equivalent to 

LTHT f = -{THT^fL & LTH = -THVLT (T^LT)H = -H\T^LT). 

s s 

Now it is easy to establish that the conditions are sufficient ("=>"): 

(1) Setting L := t N and S := T l T gives SS = T l TT^T = +l N . Thus S = S* is 
unitary, complex symmetric and satifics SH = —W S. 

(2) Setting L := J and S := T l JT gives SS = T l JTTUf = -1 N by J 2 = -l N . 
Thus S = —S is unitary, skew- symmetric and satifies SH = —H l S. 

Moreover the conditions are also necessary ("•<=") by Lemma fTTl because with 
appropriate respective unitaries T 

(1) for L = In any symmetric unitary matrix S can be written as S = T'T; 

(2) for L — J any skew-symmetric unitary matrix S can be written as S = T l JT . 

■ 

In the context of filtering simple subalgebras, Lemma 1121 can be turned into 
the powerful Algorithm 3. It boils down to checking a system of homogeneous 
linear equations for solutions S satisfying SH V = —HlS for all iH v G £ simul- 
taneously: if S is a solution with SS = +1, the subalgebra t of su(N) generated 
by the {iH„} is conjugate to a subalgebra of so(N), while in case of SS = — 1, 
t is conjugate to a subalgebra of sp(N/2). 
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Remark 13. By irrcducibility of t (via Algorithm 2), those subgroups generated 
by I C su(2") with a unitary representation equivalent to its complex conjugate 
are limited to orthogonal and symplectic ones: it follows from Schur's Lemma 
that SS = ±1 arc in fact the only types of solutions for SH = —H S with iH £ t, 
as nicely explained in Lem. 3 of Ref. |99j . Lemma [5] (of this paper) explains why 
these solutions are mutually exclusive. Due to the irreducibility, the matrix S is 
unique up to a scalar factor c £ C with cc = 1 . 

Conjugation to the symplectic algebras has also been treated in Ref. [84] by 
solving a system of linear equations, while Ref. |101j resorted to determining 
eigenvalues for discerning the unitary case from conjugate symplectic or orthog- 
onal subalgebras. — The results can be summarised and extended as follows: 

Theorem 14 (Candidate Filter I). Consider a set of Hamiltonians {iH^} 
generating the dynamic algebra { C su(N) with the promise (by Algorithm 2 
and e.g. due to a connected control system) that { C su(N) is given in an ir- 
reducible representation and t is simple. If in addition Algorithm 3 has but 
an empty set of solutions, then fi is neither conjugate to a simple subalgebra 
of Sp(N/2) nor of so(N). In particular, I is none of the following simple Lie 
algebras: C7, eg, fi, or 02. 

Proof. The cases so(N) and sp(N) are settled by Lemma IT2"1 The cases 67, t&, 
f4, and Q2 follow from the elaborate classification of Malcev [52] (see also, e.g., 
88,102.38] and Theorem H5] in Appendix IC.3[) . as an irreducible representation 
of Cs, fi, or Q2 is always conjugate to a subalgebra of so(N), while an irreducible 
representation of is conjugate either to a subalgebra of so(iV) or of sp{N/2). 



7.2. Unitary Subalgebras. It also follows from Malcev [53] (again, see also [551 
I102ll3"5] and Theorem |35J in Appendix IC.3[) that only the subalgebras su(£ + 1) 
^ 2), so(4£+2), and can have unitary representations. One can immediately 
deduce from the Tables [1] and [5] the following 

Corollary 15. The Lie algebras su(2") do not possess (proper) unitary, irre- 
ducible simple subalgebras ifn £ {1, 2, 3, 5, 7, 9, 11, 13, 15}. In these cases (n ^ 1) 
and under the conditions of Theorem \lJ\ Algorithm 3 provides a necessary and 
sufficient criterion for full controllability. M 

We checked by explicit computations that z§ does not occur as a unitary, 
irreducible simple subalgebra of su(2 ra ) for n < 100, i.e. for qubit systems with 
up to 100 qubits. Thus one might conjecture that z§ does not occur as a unitary, 
irreducible simple subalgebra for qubit systems in general. 

We present an example of a control system whose dynamic algebra is a 
(proper) unitary subalgebra of su(2 4 ): 

Example 16. Consider a bilinear control system on su(16) with four subsystems 
given by su(2) ©su(2) ©su(2) ©su(2). The local dynamic algebra is given by 

(iXIII, zYIII, zZIII, illLX, THY, TIIZ) Lic . 

In addition, we have a drift Hamiltonian H d = XXII + YYII + IXXI + IYYI + 
IIXX + IIYY (Heisenberg-XX interaction). The control system 
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is connected and acts irreducibly. The dynamic algebra t = so (10) is simple and 
a (proper) unitary subalgebra of su(16). 



7.3. System Algebras Comprising Local Actions su(2)® n . We now discuss the 
set of local unitary transformations SU(2)®" C SU(2 n ) and its Lie algebra 
su(2)® n C su(2") where both are given in their respective standard representa- 
tion, i.e. as n-fold Kronecker product and ?i-fold Kronecker sum (see Sec. 01 

5u(2)§su(2)§ • • ■ §5u(2). 

What is the classification of su(2)®" w.r.t. symplectic, orthogonal, and unitary 
subalgebras? We obtain from Thm. 3 of Ref. [TT] (see also [ITJ51IT04] ) : 

Lemma 17. For the algebra su(2)®" given in its (irreducible) standard repre- 
sentation there are two cases: (1) if n is odd, it is a symplectic subalgebra of 
su(2 n ) in the sense of being conjugate to a subalgebra of sp(2™ ), and (2) if 
n is even, it is an orthogonal subalgebra in the sense of being conjugate to a 
subalgebra o/so(2"). 

Proof. Let (j> denote an irreducible representation of a compact Lie group G. 
Then for the Frobenius-Schur indicator (Chap. IX, App. II. 2, Prop. 4 of Ref. [75] ) 
one finds 

[ — 1 (j) is a symplectic representation 

/ Tr[</> 2 (g)] dg = < +1 <=> is an orthogonal representation 
is a unitary representation 

Let ip denote the standard representation of the Lie group H = SU(2)®". Ref. [TT) 
proves that J H Tr[ip 2 (h)] dh = (-1)". ■ 

If the subsystems of a control system are independently fully controllable then 
it follows from Lemma [T71 that some cases can be excluded: 

Lemma 18. Assume that the dynamic algebra t C su(2") is irreducible and 
simple, and that the subsystems su(2) are independently fully controllable [i.e. t 3 
su(2) ©su(2) © ••• ©su(2)/. If n is odd [resp. even] then 6 is not an orthogonal 
[resp. symplectic] subalgebra. 

Proof. We remark that f) = su(2) ©su(2) © • • • ffisu(2) is given in an irreducible 
representation. It follows from the discussion prior to Lemma |H] that \) has an 
invariant (non-degenerate) skew-symmetric bilinear form and no invariant (non- 
degenerate) symmetric bilinear form if n is odd. Therefore, the dynamic algebra 
t cannot have an invariant (non-degenerate) symmetric bilinear form and the 
Lemma follows for odd n. The case of even n is similar. ■ 

Unfortunately 4 , Lemma [TH] is no longer true if the dynamic algebra contains 
only a non-zero subset of the local operations: 

4 In Thm. 3 and 4 of the conference paper Ref. 1981 we incorrectly gave more general results 
for dynamic algebras which contain a non-zero subset of the local operations. But in light of 
Examples I16l and ll9l thc more general results in Ref. 98 arc not correct, as the non-zero subset 
is in general not given in an irreducible representation. 
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Example 19. Consider a bilinear control system on su(8) with three subsystems 
given by su(2) © su(2) © su(2). The local dynamic algebra is (iXIl, iYII, iZII) Lic . 
In addition, we have a drift Hamiltonian H d = XXI + YYI + IXX + IYY. The 
control system 




is connected and acts irreducibly. The dynamic algebra t = so(7) is simple and 
an orthogonal subalgebra. We emphasise that as a consequence of Lemma [15] 
this would have not been possible if t 2 su(2) ffisu(2) ©su(2). 



7.4- A Necessary and Sufficient Symmetry Condition. In this subsection we 
present a necessary and sufficient symmetry criterion for full controllability 
of control systems contained in su(N). To this end, we introduce some addi- 
tional notation: Assume that <fi is a representation of a compact Lie algebra 
of dimension N. The tensor square <fi® 2 :— <j> ® Ijv + Ijv ® <t> decomposes as 
ifi® 2 = Alt 2 (/) © Sym 2 c/), where the alternating square Alt 2 </> and the symmetric 
square Sym 2 <^> are the restrictions of 4>® 2 to the antisymmetric and the symmet- 
ric subspace, respectively. More details on this notation is given in Appendix iDl 
We arrive at 

Theorem 20. Assume that t is a subalgebra ofsu(N) and denote by <P the stan- 
dard representation ofsu(N). Then, the following are equivalent: 

(1) t = Sll(N). 

(2) The restrictions o/Alt 2 <£ and Sym 2 <? to the subalgebra t are both irreducible. 

(3) The commutant comm$»2 (t) w.r.t. the tensor square <I>® 2 has dimension two. 

Proof. (1) =>■ (2) follows by Theorem 1551 in Appendix D. We prove (2) (1). As 
the restriction (Alt ^)|{ of Alt $ to t is irreducible, we get that the restriction 
<P\i of # to t is also irreducible. Otherwise, <P\t = <\>\ © (j>2 would be reducible 
and, as a consequence, (Alt 2 #)| t = Alt 2 Oi©0 2 ) = Alt 2 ^iffi(^i©0 2 )ffiAlt 2 2 
would also be reducible (which is impossible) . Lemma [3] implies the centraliser 
t' of t in su(N) is trivial, thus by Lemma 0] Ms semisimplc. Now (1) follows 
by Theorem 1551 in Appendix D. Moreover Thm. 1.5 of Ref. |105j says that the 
dimension of the commutant of a representation <f> is given by mf where the 
rrii are the multiplicities of the irreducible components of 4>. As we consider the 
representation (^® 2 )| f = (Alt 2 <£)| t © (Sym 2 ^)| { , the equivalence of (2) and (3) 
readily follows. ■ 

We now show that condition (3) of ThcorcmrSHlcan be easily tested using a set 
of Hamiltonians {iH„} generating the dynamic algebra t C su(N). Therefore, we 
prove that the commutant of {(iHv) <£) 1 m + 1 n ® (iH v )} is equal to comm<p®2 (t). 
Obviously, the latter commutant is contained in the former. Let s e gI(iV 2 ,C) 
be an element of the former commutant. Then by Jacobi's identity [[a, b],s\ + 
[[&, s], a] + \[s, a], b\ = 0, s commutes with all commutators 



[(iH„) g> Ijv + In ® (iH„), [iH^) <g> 1 N + t N ® (iH^)} 
= [iH u , iH^} ®t N + t N ® [iH v , iHfj], 
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and by induction, s is also contained in the latter commutant. 

Together with Theorem 1201 we thus obtain a necessary and sufficient sym- 
metry condition for full controllability as a theoretical main result: 

Corollary 21. Consider a set of Hamiltonians {iH u \ v = d, 1, 2, . . . } generat- 
ing the dynamic algebra t C su(N). The corresponding control system is fully 
controllable in the sense t = su(N), if and only if the joint commutant of 
{(iH u ) ® Ijv + 1jv <® (iH u ) \ v = d, 1, 2, . . . } has dimension two. I 

In spite of the beauty of simplicity of this result, from an algorithmic point of 
view the above symmetry condition is currently not appealing: In Corollary 1211 
one would have to compute the commutant of N 2 x N 2 matrices as compared 
to N x N matrices in the test for the lack of symmetry in Algorithm 2. Thus 
the complexity of testing for Corollary [21] would be the square of the complexity 
of Algorithm 2. Even in moderately-sized examples one has to save computer 
memory by methods of sparse matrices due to the larger matrices. In larger 
examples, testing for Corollary[5T]gets impractical. Yet compared with potential 
conditions involving even higher tensor powers, one should consider Corollary [21] 
as a fortunate incidence. 

In order to characterise the commutant of Corollary [5T] in further detail, we 
introduce an iV 2 x N 2 permutation matrix Kn^n also known as commutation 
matrix |106U107| . Let e a denote the vector such that (e a )b — S a ,b with a, b € 
{1, . . . , iV 2 }. We define Kn,n by the permutation Kn,n • e a = e^ b > where one 
has ir(N ■ i + j + 1) = j ■ N + i + 1 and i, j <G {0, . . . , N — 1}. The commutation 
matrix operates on the vec-representation |106j of an N x N matrix A as the 
transposition operator: Kn : n ■ vec (A) = vec {A*). 

Lemma 22. The commutant of ' {(iH v ) ® Ijv + 1/v ® always contains the 

elements 1^2 andK^,N- 

Proof. As the identity matrix 1^2 always commutes, we have only to prove that 
Kn,n is contained in the commutant. Sec. 3 of Ref. |107j says that Km,n{A® 
B) = {B® A)K NtN for all N x iV matrices A and B and thereby K NyN (A®B + 
B®A) = (A®B+B®A)K N)N . In particular one finds K NtN (A®t N + t N ®A) = 
(A ® Ijv + 1/v ® A)Kw y N and the Lemma is proven. ■ 

The operator has two eigenspaces (see Sec. 4.2 of Ref. [ 1 T] ) : The first 

one is given by the symmetric subspace (i.e. 'bosons') and has the eigenvalue +1 
with multiplicity N(N + 1)/2. For even TV, the permutation-symmetric subspace 
is equivalent to the Lie algebra sp(N/2). The second one is given by the anti- 
symmetric subspace (i.e. 'fermions') and has the eigenvalue —1 with multiplicity 
N(N — l)/2. The permutation-antisymmetric subspace is equivalent to the Lie 
algebra so(N). 

The methods of this subsection thus shed new light on the symplcctic and 
orthogonal subalgcbras (see Subsection 17. ip . Prop. 3.5 of Ref. |108j (see also 
p. 446 of Ref. |109j ) says that an irreducible representation of a compact 
simple Lie algebra q is either symplcctic or orthogonal if and only if its tensor 
square 4>® 2 contains the trivial representation of g exactly once. In particular, 
the irreducible representation <f> is symplectic (resp. orthogonal) if the trivial 
representation occurs exactly once in Alt <fi (resp. Sym <fi). A similar condition 
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is given by Prop. 4.2 of Ref. |108j : An irreducible representation </> of a compact 
simple Lie algebra g is either symplectic or orthogonal if and only if its tensor 
square (p® 2 contains the (irreducible) adjoint representation of a at least once. 



8. Simulability 

Simulating quantum systems [Tll5l lll0j is a promising mid-term perspective, be- 
cause the accuracy demands are easier to come by than the 'error-correction 
threshold' for actual quantum computing. Another practical advantage lies in 
the fact that sometimes the simulating systems allow for separating control pa- 
rameters in the analogue that in the original (be it classical or quantum) cannot 
be tuned independently. 

This section exploits that the dynamical algebra captures all the key proper- 
ties of the dynamical system to be studied. More precisely, the question whether 
(and to which extent) one quantum system can simulate another one can be an- 
swered by analysing the Lie-subalgebra structure of systems with a given dimen- 
sion. Recently Kraus et al. have explored whether target quantum systems can be 
universally simulated on translationally invariant lattices of bosonic, fermionic, 
and spin systems |17| . Based on the branching diagrams of simple subalgcbras to 
su(N), here we take a more general approach pursuing the question which type 
of quantum system can simulate a given one with least overhead in state-space 
dimension. In particular, we also allow for effective many-body interactions to 
be simulated by pair-interactions. — To this end, the reader may wish to resort 
to the more general notion of tensor-product structures in Appendix [Al first . 

In quantum simulation, one of the first natural questions to ask is whether 
and under which conditions a controlled quantum dynamical system S a can sim- 
ulate another (controlled or uncontrolled) dynamical system Sb given as bilinear 
control systems with p = a, b on density matrices p M 

m 

Pfl (t) = -i [(ftf + J2 u$(t)H?) , Pli (t)] with ^(0) = Pllo . (15) 

The dynamic Lie algebras t a and 6f, are given by the respective Lie closures as 
fi p := (iHg,iH?\j = 1,2,..., m) U e (16) 

thus entailing the reachable sets take the form of K„-subgroup orbits as in 
Eqn. Q) 

react) (pao) := {K a p ao K\ \ K a G K a := cxp« a } and (17) 
teacfj(/5 6o ) := {K h p ho K\ \ K b G K b := expfi,,} . (18) 

An obvious requirement is that for any initial state pbo of system Sf, leading 
to the dynamics Pb(t) £ teacf)(/9& ) there is an initial state p ao of system S a such 
that under the dynamics of S a one has 

Pb(t) G ttach( Pa o) Vt>0. (19) 



This requirement is obviously fulfilled by the following sufficient condition: 
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Proposition 23. A dynamic bilinear control system E a with dynamical algebra 
t a can simulate another dynamic system Eb with dynamical algebra tb ift a 2 fy>- 

Proof. Clearly t a D tb implies K a D Kb and thus teact)(p ao ) D reacf)(/jf, ), which 
in turn ensures that Eqn. (|19p is fulfilled for any choice of initial states. ■ 

In particular, if system Eb is uncontrolled it can be simulated if its drift Hamil- 
tonian Hq can be simulated, i.e. provided iH^ £ 6„. 

Two dynamic bilinear control systems E a and Eb are said to be dynamically 
equivalent independent of the respective initial states p^o if and only if they can 
mutually simulate one another, i.e. if t a 3 h and tb 3 £ a so t a = h (up to 
isomorphism) . 

Remark 24- It is important to note that in the special case of pure states, where 
by construction p(t) = p 2 (t), it suffices that, e.g., a system E a has the dynamic 
Lie algebra 6 a = sp(N/2) in order to simulate system Eb with %b = su(N), be- 
cause the unitary orbit of any pure state po = \ip) (ip\ coincides with its symplectic 
orbit for N even 

O SV (N)(\'>P)(->P\) = SpW2) (\ip}(iP\) y\^)£H. (20) 

This is equivalent to a well-known result stating that for N even, a system is pure- 
state controllable as soon as its system algebra encapsulates the symplectic one 
|36j . — Since we are interested in general results beyond pure states, the notion of 
full controllability maintained in this work is full operator controllability unless 
specified otherwise. Also for simulability we do not confine the state space to 
pure states henceforth. 

Proposition 25. Consider two dynamic systems E a and Eb whose respective 
dynamic Lie algebras t a and tb shall be irreducible over a given Hilbert space 7-1. 
Then E a simulates Eb irreducibly and with least overhead in the very 1-L given, 
if any interlacing system Ei with irreducible algebra 6, fulfilling 

ta 2 ti 3 h (21) 

enforces ( up to isomorphism ) 6j = { a or ti = tb or trivially both. 

Caveat. Note that the term 'with least overhead' crucially depends on the Hilbert 
space given a priori: Thus there may be extreme realisations. For instance, in a 
fully controllable system of say 14 qubits with dynamic algebra su(16 384) there 
is an irreducible way to simulate a fully controllable su(4)-system of two qubits 
(or just a single spin-| with control over all multipole moments) with 'least 
overhead' in su(2 14 ), see the penultimate entry in Tab. [2] Realisations of this 
type may not be very useful in practice, yet relate to the context of code spaces. 

Here, we have dealt with quantum simulation of unobserved control systems. 
Now we illustrate the above findings by examples. Later, in Sec. 1121 we will give 
an outlook on a weaker notion of quantum simulation of observed control systems 
with respect to expectation values by given sets of observables. 
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Fig. 2. Heiscnberg-XX spin chains with n spins- ^ and odd-order orthogonal system algebras 
so(2n + l) require one locally controllable qubit at the end. A full series can be constructed, the 
first examples being so(5) — 5(3(4/2), so(7), so(9), and so(ll). For n = 1 one gets so(3) = su(2). 



9. Worked Examples 



9.1. Dynamic Systems with Orthogonal Algebras. Take spin chains of n spins-j 
with Heisenberg-XX (and XY) interactions and a single locally controllable qubit 
at one end. These instances serve as convenient topologies to simulate a full series 
of odd-order orthogonal algebras so(2n+ 1) for n qubits. The first instances are 
shown in Fig. [5J 

Proposition 26. Heisenberg-XX chains of n spin-^ qubits (n>l) and a single 
locally controllable qubit at one end give rise to the dynamic system algebras 
So(2n + 1) as irreducible subalgebras embedded msu(2"). 



Proof. In view of later applications, the proof is kept constructive. For better 
readability, let x, y, and z denote Pauli matrices. 

First, as a foundation for induction, the case n = 2 can be settled by direct 
calculation to verify 



i{xl, yl, (xx + yy))u c = i{xl, yl, zl, xx, yy, xy, yx, zx, zy, lz} = p so(5) , (22) 



where the final identity can be corroborated by Algorithm 3 as will be illus- 
trated in Eqn. (|2"3"j) below. 

Second, for the induction from (n — 1) to n, where the drift Hamiltonian 
is extended by the final Hciscnbcrg coupling between the qubit pair (n — l),n 

to take the form Hq := X)l-=i ^fc^fc+i + VkUk+i; observe that all the algebra 
elements for n — 1 qubits re-occur. Upon twice commuting with zl ■ ■ ■ 1 arising 
at the controlled end, the first pair coupling term X1X2 + 1/12/2 can be recovered: 

a< ii Zl (iJ2kZl x kXk+i+ykVk+i) = adi 2l (-i(2/i.T2-xiy 2 )) = -i {x\x 2 +2/12/2) and 
then by virtue of Eqn. (|2"2"T) also lzl ■ ■ ■ 1 and thus recursively all the terms in 
the Lie closure at the stage n — 1. 

Third, once having embedded the (n— l)-qubit algebra into the n-qubit sys- 
tem, the induction boils down to including the coupling term x n -\x n + y n -\y n , 
which takes 1 • • ■ lz n _il to 1 ■ ■ • lz n . Writing braces { y } whenever one has the 
choices {x,y} one gets the following complete list: 
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Fig. 3. n-Spin-g Heisenberg-XX chains with even-order orthogonal system algebras so(2n + 2) 
result by allowing just two locally controllable qubits at the ends. A full series can be con- 
structed, the first examples of which are shown so(6) = su(4), so(8), so(10), and so(12). For 
n = 2 one gets so(6) = su(4) as a fully controllable two-qubit system. 



The Pauli-basis elements for so(2n + 1) 



2 terms 
n terms 




{£}!•••! 
zl---l 


lzl • • • 1 etc 




2 terms 
4(n-l) 


terms 






■ ■ 1 etc 


2 terms 
4(n - k 


■f 1) terms 


il} z --- z {y} k 1 --- 1 


l{y} 2 Z---Z 


{JWi 1 '" 1 * 


2 terms 
4 terms 




zz --- z iy} n 
{y} z --- z { X y} n 







Finally counting terms gives a total of 2n+n + 4 ^l = \{n — k) = 3n+4 X)fe=i ^ = 
3n + 2n(n — 1) = 2?i 2 + 77 = dimso(2n+ 1) = n(2n+ 1) elements to span the basis 
of the Hamiltonians H u generating (iH v )i& e = So(2n + 1). So for all n-spin-i 
Hciscnbcrg-XX chains controlled locally at one end we have obtained a construc- 
tive scheme to determine irreducible representations of their respective dynamic 
Lie algebras so(2n + 1) in terms of Pauli bases. ■ 



In contrast, n-spin-i chains with Hciscnbcrg-XX interactions and two inde- 
pendently controllable qubits, one at each end, provide a realisation of a series 
of even-oder orthogonal algebras so(2n + 2) for n qubits, the first examples being 
shown in Fig. [31 



Proposition 27. Heisenberg-XX chains of n spin-^ qubits (n > 2 ) and two 

individually locally controllable qubits, one at each end, give rise to the dynamic 
system algebras so(2n + 2) as irreducible subalgebras embedded in su(2"). 



Proof. The constructive proof follows in entire analogy to the one of Proposi- 
tion however, the local controls at the second end imply that the Lie closure 
comprises each term occuring in the above list also read from right to left thus 
duplicating the first line in each category from two terms to four terms. Since 
the second lines in each category already comprise the reverse terms, one obtains 
the following complete list of elements: 
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The Pauli-basis elements for so(2n + 2) 


4 terms 


{?}!•••! 


!•••!{?} 


n terms 


zl • • • 1 


Izl • • • 1 etc 


4 terms 




\- ■■!{%} z 


4(n — 1) terms 


{vHi/}*"' 1 


1 { y } { y } 1 ' ' ' 1 etc 


4 terms 




1 --- 1 { X y} n -k z --- z 


4(n — k + 1) terms 




■1 l{y} 2 z---z{^} fc+1 l---l etc 


4 terms 


z --- z { X y)n 




4 terms 


{ X y} z --- Z {'t}n 




1 term 


zz ■ ■ ■ zz 





Finally, by the commutator [(2 



= {S}*i---i),(i---i{?} f 



5)] with 



fc = n—k', the longitudinal spin-order term Z1Z2 ■ 



terms, one arrives at a total of 4n + n + 1 + 4 Y^j=\ ( n — i) = 5n + 1 + 4 
57i + 1 + 27i(ti — 1) = 2n 2 + 3n + 1 = dim5o(2n + 2) elements. Thus also for all 
n-spin-i Hciscnberg-XX chains individually controlled locally at the two ends we 
have provided a constructive scheme to determine irreducible representations of 
their respective dynamic Lie algebras so(2n + 2) in terms of Pauli bases. ■ 

In both instances of Heiscnbcrg-XX chains controlled locally at one end [Fig. [5] 
with so(2ti+1)] or at two ends [Fig.[3]with so(2n+2)] there are convenient Cartan 
decompositions q = £©p: the f-parts consist of per- antisymmetric matrices, while 
the p-parts comprise the per- symmetric matrices, recalling that per-symmetry 
relates to reflection at the minor diagonal. In both of the above listings, the 
respective subalgebras J to so(2n + 1) or so(2n + 2) encompass the Hamiltoni- 
ans with odd numbers of z-terms, while the respective subspaccs p contain the 
elements with even numbers of z-terms (including zero z-terms). 

For illustration, in the first example, i.e. the two-qubit Hciscnberg-XX chain 
of Fig. [5J the transformation matrix S satisfying SH + = according to 
Algorithm 3 is given by 



listed last arises. Counting 
1 . 

1 J = 



Si 







-1 



00+1 
+1 
-100 
000 



■h 



(23) 



Here S\Si = —1 reconfirms so (5) = sp(4/2). 

As a second example, for both of the three-qubit cases in Fig. [5] and Fig. [3] 
corresponding to so (7) and so (8), Algorithm 3 provides 



/0000000+l\ 
000000 -1 \ 




-1 
0-10 

















00+100000 



(I 



1 



Vu --l u w w w 
+ 1 








(24) 
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(a) (b) 
® — ® — ® ® — ® — © — ® — ® 
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Fig. 4. Quantum systems with dynamic Lie algebras sp(4) [see (a)] and sp(16) [see (b)] as 
examples of a series of Ising chains of n = 2k + 1 qubits with positive ZZ coupling terms on one 
branch and negative couplings on the other. They give rise to the dynamic algebras sp(2 n— 1 ) 
irrcducibly embedded in su(2 n ), respectively. The limiting case k = gives sp(l) = su(2) as a 
single fully controllable qubit. 



where 5*25*2 = +1 shows the orthogonal type of the respective irreducible repre- 
sentations. 



9.2. Dynamic Systems with Symplectic Algebras. Based on the smallest exam- 
ples of qubit systems with Ising-ZZ interactions shown in Fig. 01 even on the basis 
of collective controls one may construct a full sequence of n spin-^ chains with 
n odd, the dynamic system algebras of which are the symplectic ones sp(2 n_1 ). 
Note again that the bilinear control systems with symplectic system algebras are 
pure-state controllable [36], whereas they fail to be fully operator controllable. 

Proposition 28. Ising-ZZ chains ofn = 2k+l spin-^ qubits (k > 1) including k 
pairs of qubits which can be controlled simultaneously and one qubit in the middle 
of the chain which can be controlled independently as in the first row of Fig.^give 
rise to the dynamic system algebras Sp(2™ _1 ) = sp(2 2fc ) as irreducible subalgebras 
embedded in su(2 n ) — su(2 2fe+1 ). We obtain the same dynamic algebras when all 
qubits can only be controlled simultaneously as in the second row of Fig. [^J 

Proof. We focus on the dynamic algebra tk corresponding to the case when all 
2fc+l qubits can only be controlled simultaneously as in the second row of Fig.UJ 
We denote by tk the dynamic algebra corresponding to the first row of Fig. |U 
We use the notation 

Xj :=I-IXI-I, Yj :=I - IYI - I, and Zj := I • • • I Z I ■ ■ • I (25) 

j-l n-j j'-l n-j j-1 n-j 

to denote the operators which act, respectively, as X, Y, and Z on the j-th 
qubit and as the identity on all other qubits. We remark that the statements 
of the Theorem can be directly verified for k € {0,1}. We organize the proof 
in steps: first we prove that lk-i ^ tfc, second we prove that = later we 
show that £fc is given in an irreducible (third step) and symplectic (fourth step) 
representation, and in the end we prove that ik is not a proper subalgebra of 
sp(2 n_1 ) = sp(2 2fe ). Recall, that is generated by the operators 

2fe+l 2fc+l / fe 2k \ 

A = "I E X i' h = ~5 X! Y J' h = ~i I ^2 Z i Z 3+ 1 ~ Z j Z i+i • 

J=l 3=1 y'=l 3=k+l / 
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The corresponding algebra tk-i on 2fc — 1 qubits can be embedded into 2k + 1 
qubits using the operators 

2k 2k ( k 2k-l 

•9i = -| 5Z X J' - 92 = ~f X Y ;m 93 = _ 2 X] ZjZj+i - X] Z J' Z i+ 1 

J=2 j=2 \j=2 j=fc+l 



We compute repeated commutators of / 3 with /i. In the first two iterations, we 
get h = [fa, fi] = -§(Ej=i [YiZj+i + ZjY i+1 ] - £™ fe+1 [Y;Z j+1 + Z,Y J+1 ]) 
and / 5 = [/ 3 ,/ 4 ] = -fhXj - 2Ej= 2 Xj - X 2fe+1 - 2(£ } =2 Z^X^+i - 
ZfcXfc + iZfc + 2 + Ej=fc+ 2 Zj-iXjZj+i)]. Repeating this process, we obtain the ele- 
ment.f 6 = [/ 3 ,/ 5 ] = -| [-Y 1 Z 2 +Z 2fc Y 2fc+1 -4(E J =2 Y,Z j+1 -E ^ fe+1 Y J -Z j+1 + 
E.?=i ZjY J+ i — E^=k+i ZjY 3+ i)]. Finally, we compute the next element ,/V = 
[/3j /e] = — |(Xi + 8 Ej=2 Xj + X 2 fc + i + 8 Ej= 2 Zj-iXj-Zj- + i — 8ZfeX/ £+ iZ/ £+2 + 
8E?=k+ 2 Z,--iX,-Z i+1 ). We obtain that /„ = -(4/ 5 + / 7 )/3 = -^Xj + X 2fc+1 ) 
and 5i = /i - / 8 . The proof for f 9 = -|(Yi + Y 2fc+ i) and g 2 = / 2 - / 9 is 
similar. We compute a few more commutators: First, we set /io = [/3,<7i] = 
-i(Ej= 2 [YjZj+i + ZjYj+i] -E -=fcii [Y,Z J+1 + Z J Y J - +1 ] + Z 1 Y 2 -Y 2fc Z 2fe+1 ). 
The other commutators are f n = [/g,/io] = -§(-YiY 2 + Y 2fc Y 2fc+ i), f 12 = 
[/i,/n] = -|(-ZiY 2 - Y!Z 2 + Z 2fc Y 2fc+1 + Y 2fc Z 2fc+1 ), and / 13 = [A,/i 2 ] = 
- 2 -(2Y 1 Y 2 -2Z 1 Z 2 -2Y 2fc Y 2fc+1 +2Z 2fc Z 2fc+1 ). It follows that f u = -|/i 3 -/n = 
-|(Z X Z 2 - Z 2fc Z 2fe+ i) and g 3 = f 3 - f u . We obtain t fc _i C completing the 
first step of the proof. Relying on the form of f% and fg we can prove by induc- 
tion that tk = tk (second step). Assuming by induction that tk-i is irreducibly 
embedded on 2k — 1 qubits, we obtain that the centralizer of tk-i (embedded 
on 2k + 1 qubits) is given by all operators O which operate only on the two 
outer qubits. But the generators /i, / 2 , and / 3 of tk do not simultaneously com- 
mute with operators O. Therefore, tfc is irreducibly embedded on 2k + 1 qubits 
(third step) . We switch to a new basis by reordering the qubits according to the 
numbers in the figure: 

^ +ZZ +ZZ -zz -zz 

••• © — © — © — © — © ••• 



In this basis, we can provide a matrix 



/0 1\ 
0-100 
0-10 

\l 0/ 



which satisfies SH + = for all elements iH of tk- In particular, we have 
SS = — 1 2 ». This can be readily verified on the generators for k € {0, 1}. Using 
our commutator computations we obtain that tfc = / 2 , fa}hie is equal to 
(<7i;S2, <?3 7 /s, /9j /i4)Lie- Thus we can prove 5i7 + H*S = by induction on fc: 
Assuming the equation holds for <7i,<7 2 , and g 3 (i.e. for k — 1), we need to prove 
that it also holds for /s,/g, and /14 which are respectively given in the new 
basis by - %(X 2k + X 2fc +i), -|(Y 2fc + Y 2fc+ i), and -§(Z 2fe _ 2 Z 2fc - Z 2fe _iZ 2fe+ i). 
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(a) 



(b) 



+ZZ +ZZ -zz -zz 
® (D © ® (A) 



(® 7®~\ 



_ +ZZ +ZZ -zz -zz 
(A) ® ® ® (A) 

_ +ZZ +ZZ -zz -zz 

® — ® — ® — ® — ® 



( ® ~ ~ " w " I _© / ~ ®~ " _ _ 



Fig. 5. Quantum systems with dynamic Lie algebra sp(16). (a) Examples with pairwise Ising- 
ZZ couplings and (b) examples with three-body ZZZ-intcractions. 



But this can be directly checked on the four outer qubits using 5*2 = M ® M. 
As 6ft is given in an irreducible representation, the matrix S is unique up to a 
scalar factor. This shows that 6^ is given in a symplectic representation and that 
tfc Q sp(2 2k ) (fourth step). Staying in our new basis, we prove that tk contains 
the elements Pj := -f(XjZ j+1 - Zj-X^+i) and Q 3 := -| (X 3 Y j+1 - Y 3 -Xj+i) 
for all even j € {2, ...,2k} by induction on j. This can be readily verified 
for j = 2 considering ti C t k . Assuming that tk contains the elements Pj-2 
and Qj-2 for j < k, we show that it also contains the elements Pj and Qj. 
Recall that tk contains the elements f%, /g, and /14. In addition, the elements 
vi = -f (X 2fc - 2 +X 2 fe_i) and v 2 = -%{Y 2 k-2 + Y 2 fc_i) are contained in fi fe . But 
one can directly check on the four outer qubits that Pj and Qj are contained 
in the algebra m = (/g, fa, fu, Pj-2, Qj-2, vi, v 2 )ue = so(2 4 ). Assuming that 
lj = sp(2 2 ^') holds for j < k, it follows that sp(2 2fe " 4 )®so(2 4 ) C t k C sp(2 2k ). As 
sp(2 2fc - 4 )(giso(2 4 ) is a maximal subalgcbra of sp(2 2fc ) (see Thm. 1.4 of Ref. [88] ) 
and € tk is not of product form, we obtain by induction that { = sp(2 2fc ). ■ 

Note the Cartan decomposition in the antisymmetric Ising chains of Fig.0]can 
be taken with respect to the joint permutation of the qubits in the two branches 
with positive and negative ZZ couplings: the t-part consists of all terms with odd 
numbers of Pauli operators deviating from the identity, while the p-part collects 
the ones with even numbers. 

As a third example, consider the first Ising chain in Fig. @] corresponding to 
sp(8/2). Here Algorithm 3 gives 



/ 

/ n 



s. 



-1 
000 +1 0000 
00000+100 
0-1000000 
000000 +1 



\ 



0-10 





10 



(26) 



u u u u — 1 u u u 

\+l 0000000/ 



with S3S3 = — t underscoring the irreducible representation is symplectic. 

Moreover since all the dynamic systems in Figures [5j a ) ari d [5][b) share the 
same system algebra sp(16), any two can mutually simulate eachother by Propo- 
sition 1231 So remarkably enough, the spin chains in Fig. [SJa) can simulate the 
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Fig. 6. n-Spin-j Heisenberg-XX chains with n {b3} and only one locally controllable 
qubit at the second position have orthogonal system algebras so(2") if (n mod 4) G {0, 1} and 
symplectic system algebras sp(2" — 1 ) otherwise. A full series can be constructed for n > 3, and 
the examples shown for n £ {4, 5, 6, 2} correspond to so(16), so(32), sp(32), so(5) = sp(4/2). In 
the single case of n = 3, central symmetry arises, which makes the respective algebra reducible. 



effective threc-qubit ZZZ-interactions shown in Fig.[5jb). In particular, note the 
lowest instance in Fig.jSJa): even only the collective local controls on all the qubits 
suffice to generate the three-body interactions with full local control shown at 
the top of Fig.(SJb). In turn, it may be astonishing at first sight that the system 
on top of Fig. [Hh) does not provide more dynamic degrees of freedom than the 
collective system at the bottom of Fig. 03a), where the simulating power roots 
in the opposite signs of the couplings. 



9.3. Dynamic Systems with Alternating Orthogonal and Symplectic Algebras. 
Based on the smallest examples of Hciscnbcrg-XX chains with one single local 
control on the second qubit as shown in Fig. El one may construct a full sequence 
of n spin-i chains, whose dynamic system algebras are orthogonal or symplectic 
depending on the value of n ^ {1, 3}. Again, observe symplectic system algebras 
ensure pure-state controllability |36| without full operator controllability. 

Quite remarkably, full local control on a single qubit suffices to get a dynamic 
algebra, where the number of dynamic degrees of freedom scales exponentially 
with number of qubits, a finding described only for full isotropic Hciscnberg-XXX 
coupling up to now [41] . More precisely, one arrives at the following: 

Proposition 29. Heisenberg-XX chains ofn ^ {1, 3} spin-^ qubits with only one 
locally controllable qubit at the second position give rise to the dynamic algebras 




so(2 n ) if (n mod 4) g {0, 1}, 
sp(2"- 1 ) if (n mod 4) € {2, 3} 



which are irreducibly embedded in su(2' 1 ). 

Proof. In the notation of Eqn. ([23)) the generators of the dynamic algebra t n can 
be written as h = -|X 2 , f 2 = -|Y 2 , and f 3 = -|(£"=i X^+x + Y,Y J+1 ). 
We remark that the statements of the Theorem can be directly verified for 
n 6 {2, 4, 5}. Assuming n > 6 from now on, we complete the proof by induction. 
We organize the proof in steps: first we prove that t n D 6 2 © ?n-2, second we show 
that t n is given in an irreducible representation, and in the end we prove that t ra 
is equal to so(2") or sp(2™~ 1 ). By computing sums of commutators we identify 
certain elements of £ n . The first elements arc f± := [^3, fi]} + [/ 2 , [/3, / 2 ]] = 
-|(X 1 X 2 +Y 1 Y 2 + X 2 X 3 +Y 2 Y 3 ), h ■= [[/ 2 ,/i],[/i,[/ 2 ,/ 4 ]]] = -|(X 1 X 2 + 
X 2 X 3 ), and f 6 := / 2 ], [/ 2 , [A ,/ 4 ]]] = — |(YiY 2 + Y 2 Y 3 ). Next we compute 
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Fig. 7. n-Spin-^ Heisenberg-XX chains with n > 2. where the first two qubits can be inde- 
pendently, locally controlled have fully-controllable system algebras su(2 n ). A full series can be 
constructed: the first examples shown correspond to the algebras so(6) = su(4), su(8), su(16), 
and su(32). 



the elements f r := [/ 4 , [/ 3 , U}} + [h, [h, [h, [/a, /*]]]] + [/i, [/s, [/i, [/a, / 4 ]]]] = 
-i(X 3 X 4 + Y 3 Y 4 ) and f a := [[/ 2l /x], [[/i,/ 2 ], [/ 3 , [/ 4 ,/ 3 ]]]] - fr + [/e, [/rja]] + 
[/ 5 , [/ 7l /a]] + [/2, [/e, [/a, [/a, /*]]]] + [/i, [h,[fi, [/a, / 4 ]]]] = -§(X 2 X 3 + Y 2 Y 3 ) 
leading to the elements f 9 := f 4 -f s = -§(X:X 2 + YiY 2 ), /io : = h~ U = 
~ I (£"=3* XjX.+i + Yj ■ Yj+i ) , and / n := / 4 + / 7 = - f (X 4 X 2 + Y 4 Y 2 + X 2 X 3 + 
Y 2 Y 3 + X 3 X 4 + Y 3 Y 4 ). By explicit computations on the first four qubits one 
can show that the elements / 42 = — |X 4 and /i 3 = — |Y 4 are contained in 
64 = (fi,f2,fu)ue Q fin- We obtain that 6 2 = / 2 , f 9 ) U e Q K and t„_ 2 = 
(/12, /i 3 , Ao>Lic C £„. Therefore, fi„ = / 2 , / 9 , /12, /i3, Ao, /s>Li e 2 *2 ffi*«- 2 - 
This completes the first step of the proof. By induction t 2 and £ rl _ 2 are given 
in an irreducible representation. Therefore, this holds also for t 2 ®t n - 2 and B n , 
which completes the second step of the proof. Using the matrices 



S 2 := 



000+1 
+1 I 

-10 



3 • — 



/ 
/ n 



+1 



+1 












0-10 









'1 
-1 0000000 



0+10 



0+10000 




\ 



(S 2 )® k and S 2k+1 = (S 2 )®( k -V 

2k+i. Relying oj 

computations in the case of j £ {2, 4, 5}, one can verify that SjH + H l S 



we define the matrices S 2k 
that S 2k S 2k = (-l) fc l 2 2 fc and S 2k +iS 2k+1 



S3. We obtain 
(— l) fc l 2 2fc+i . Relying on direct 



holds for all elements iH of tj. Assuming by induction that SjH 





H*Sj = 

holds for all elements iH of tj where j £ {2, n— 2}, we show that SnH+H 1 S n = 
holds also for all elements iH of the algebra t„ = (fi,f 2 , fg, f 12 , f 13 , fw, fa) hie 2 
fi 2 ©fi«-2 by directly verifying S^fs + fs^4 = on the first four qubits. In 
summary, we proved that t 2 © t n - 2 Ct„C so(2") or sp(2 n_1 ) depending on the 
value of n. But Thm. 1.4 of Ref. [55J says that t 2 © t n - 2 is a maximal subalgebra 
of so(2 n ) or sp(2™" 1 ). Thus, t n is equal to so(2") if (n mod 4) G {0, 1} and equal 
to sp(2 n_1 ) otherwise. This completes the last step of the proof. ■ 



9.4. Dynamic Systems with Unitary Algebras. We close the series of worked ex- 
amples by considering n-spin-^ Heisenberg-XX chains with n > 2, where the 
first two qubits can be independently, locally controlled (see Fig. [7J). This case 
was recently studied in Refs. [4Tll4l^[TTl"] . We show that these systems are fully 
controllable for arbitrary n > 2. 
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Corollary 30. Assume that the first two qubits of a Heisenberg-XX chain of 
n spin-^ qubits with n > 2 can be independently, locally controlled. Then, the 
dynamic algebras is t n =su(2 n ). 

Proof. The Theorem can be directly verified for n E {2, 3, 4, 5}. Building on the 
proof of Proposition 1291 we prove the Theorem for n > 6 by induction. We first 
show that t n 3 ?2©tn-2- From the proof of Proposition [25] it is only left to 
show that the elements X3 and Y3 can be generated. But this can be directly 
verified by computations on the first four qubits. Thus we proved that su(2") D 
fin 2 su(2 2 )§su(2™- 2 ). Thm. 1.3 of Ref. [55] says that su(2 2 ) ®su(2™- 2 ) is a 
maximal subalgcbra of su(2"). The Theorem follows immediately. ■ 

10. Fermionic Quantum Systems 

Fermionic d-level systems with any kind of quadratic (pair-interaction) Hamilto- 
nians give rise to dynamic system Lie algebras limited to subalgebras like so (2d) 
or so(2d + 1). By making use of the Jordan- Wigner transformation, which links 
the number of levels d with the number of qubits n, we show how these systems 
can be simulated by n-spin-i chains with partial local control. — For keeping 
the relation to mathematical literature, references are more extensive in this 
section. 

10.1. Quadratic Hamiltonians. To fix notations, consider the fermionic creation 
and annihilation operators f£ and f v which operate on a finite-dimensional 
quantum system of d levels and satisfy the anticommutation relations (with 
1 < P, <7 < d and {a, b}+ := ab + 6a and the Kronecker function S p _ q ) 

{fi,f q }+ = 5 p , q and {f p ,fl} + =0 = {f p J q } + . (27) 

For the p-th level of the system, /t and f p change the occupation numbers n p 
labelling the respective states \n p ) such as to give /J|0) = |1) = (J) = \\) and 
/ p |l) = |0) = (?) = ||), where by the Pauli principle n p e {0,1}, (/ p t) 2 = 0, 
and (/p) 2 = 0. Now the properties of the usual quadratic Hamiltonians (see, e.g., 

[naimiiniEiins]) 

d 

H .= Apgfpfq + Bpqfpf^+Cpqf^fq+Dpqf^f^. (28) 

p,g=l 

can be discussed in terms of their pair-interaction coupling coefficients A pq , B pq , 
C P q, and D pq seen as (possibly complex) entries of the dxd- matrices A, B, C, 
and D. Hcrmiticity of H requires A = D\ B = B\ and C = C\ while in 
addition, the commutator relations of Eqn. (|27[) imply 

d d 

H = ^ -A p qfqf p - B> vq f\fl + ^ ~^pqfqfp ~ Cpqfqfp 
p,q=l p ,q=l, p ^q 
d 

+ Y, B„(l - flf P ) + C pp (l - f p fl) , 
p=i 
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which upon identification with Eqn. (J2HJ) enforces A = —A*, B = — C*, and 
D = —D l . Finally keeping a widely- used custom (see, e.g., p. 452 of Ref. |112) 
or p. 173 of Ref. |113j ) we also assume the entries of A, B, C, and D are real. 
Summing up, A is real skew-symmetric following A = A = —A 1 = —D and B 
is real symmetric with B = B = B 1 = — C . So H of Eqn. ([38"|) can be given in 
'symmetrised' form (see, e.g., p. 2 of Ref. \11Q\ ) as 

d 

p,q=l 



10.2. Jordan- Wigner Transformation. For simplicity, first recall the map from 
the non-compact, (real) special linear algebra sI(2,R) to the compact, special 
unitary algebra su(2). The generators of s[(2,R) are given by E = (q J) = 
\ (X + i Y) , H = ( J ) = Z , and F = ( ° g ) = ± (X - iY) and the generators of 
su(2) can be chosen as iX, iZ, and iY where X=(5o)=F + E, Z = (q — i) = 
H, and Y = (9 ' Q l ) = i(F — E). Obviously, this also defines a map between 
the Lie algebras (see, e.g., p. 127 of Ref. [15] or Chap. IX, Sec. 3.6 of Ref. [76"] ) 
readily serving as a prototype for maps between non-compact normal real forms 
and compact real forms of Lie algebras (cp. [77]). 

Likewise, the Jordan- Wigner transformation [117] maps the fermionic operators 
(for 1 < p < d) 

fp = \{c P + ic p+d ) and f p = | (c p — icp+d) to the operators (30) 
c P ■= f P + fl and c p+d := i(f p - /+) . (31) 

Now c p and c p +d can be given explicitly as d-qubit operators 

c p = Z_^_ZXL^_I and c p+d = Z - ZYI - I . (32) 

p— 1 d—p p—1 d—p 

We refer to Chap. VIII, Sec. 3 of Ref. [TT8], Sec. 9.6 of Ref. PH], or Sec. 44 
of Ref. [120] . where more information on this construction can be found and 
where connections to Clifford algebras are discussed. In the context of Clifford 
algebras this construction is sometimes named after Brauer and Weyl [121] (see, 
e.g., p. 183 of Ref. [123] L 



10.3. Quadratic Hamiltonians in Qubit Form. Now one can readily apply the 
Jordan- Wigner transformation to fermionic quadratic Hamiltonians. Assuming 
that the number of levels is d, the Hamiltonian of Eqn. ([39"]) is mapped to (see, 
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e.g., Thm. VI.I of Ref. [TT5] ) 



d 



H = - Y b pp 1 • • • 1 z 1 • • • 1 



(33a) 



P =1 p-1 d-p 



(1 




(33b) 



P,9=1,P>9 <j-l p-g-1 d-p q-1 p-g-1 d-p 



(1 




) 



(33c) 



p,q=l,p>q g _l p-g-1 d-p q-1 p-«-l d-p 



This determines the dynamic algebra of a general fermionic Hamiltonian con- 
taining quadratic terms: 

Theorem 31. Let the entries of the real antisymmetric matrix A and the real 
symmetric matrix B denote the control functions of the Hamiltonian given in 
Eqn. (|33p . We assume d > 2. The dynamic algebra so (2d) of the corresponding 
control system is embedded in su(2 d ). The centraliser of the dynamic algebra is 
one- dimensional and is given by the d-qubit operator — |Z • • • Z. The embedding of 
so(2d) into su(2 d ) splits into two irreducible representations of equal dimension. 

Proof. Let td denote the dynamic algebra of the control system. The generators 
■ ■ IXZ • • ■ ZXI • • • I and -|l ■ ■ • IYZ • • • ZYI - I arise from linear combina- 
tions of Eqns. (|33bp - (|33c[) . and computing commutators with the generators 
— |l- - ■ IZI- • -I from Eqn. (|33a| reveals the generators -| I • • ■ IXZ • • ■ ZYI - I 
and — |I • • ■ IYZ • • ■ ZXI ■ • ■ I. By comparing with the (independent) proof of The- 
oreml33lit follows that td is a subalgebra of so(2d+l). The statements of the The- 
orem can be directly verified for d G {2, 3, 4, 5}. We assume by induction that the 
(d— l)-qubit operator a = — |Z • • • Z is the only element in the centraliser of td-i- 
Considering a as an d-qubit operator operating on the first d — 1 qubits we obtain 
that the centraliser of td can only contain linear combinations of elements from 
the set — |{I ■ ■ ■ IX, I • • ■ IY, I ■ ■ • IZ, Z • • ■ ZI, Z • ■ ■ ZX, Z ■ ■ • ZY, Z • ■ ■ ZZ}. The sec- 
ond statement of the Theorem follows from the fact that only the last element in 
the set commutes with all generators. We obtain from the structure of the gener- 
ators that so(2d — 4)©so(4) C t d - We remark that so(2d— 4)©so(4) is a maximal 
subalgebra of so (2d) and that so (2d) is a maximal subalgebra of so(2d+ 1) (see, 
e.g., p. 219 of Ref. [H3] or Table 12 on p. 150 of Ref. [121]). Therefore, t d is 
equal to so(2d) or so(2d + 1). But td ^ so(2d + 1) as the corresponding em- 
bedding would be irreducible, and the first statement of the Theorem follows. 
We already showed that the centraliser is one-dimensional which is equivalent to 
the fact that the commutant is two-dimensional. Theorem 1.5 of Ref. [105] says 
that the dimension of the commutant of a representation is given by mf 
where the m, are the multiplicities of the irreducible components of <f>. Thus, 
the embedding of so (2d) into su(2 rf ) splits into two irreducible representations. 
The third statement of the Theorem follows now by proving that the simulta- 
neous eigenvalues of the Hamiltonians corresponding to all generators are given 
by ±1 occurring each with multiplicity 2 d ~ 1 . This can be directly verified for 
d G {2,3,4,5}. Assuming the statement by induction for all td' with d' < d 
we obtain the simultaneous eigenvalues of Hamiltonians corresponding to the 
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algebras td-2 and £2 acting on the first d — 2 qubits and the last two qubits, 
respectively. As the eigenvalues of a tensor product of two matrices are given by 
the product of the eigenvalues of each matrix, we can prove the statement by 
induction. ■ 

The first statement of Theorem [21] is related to the fact that the canonical 
transformations of fermionic systems are given by orthogonal transformations 
o(2d) (see, e.g., p. 118 of Ref. |113j and p. 39 of Ref. |114j ). In more mathematical 
literature, the first statement of Thcorcm[2Dcan be found in Sec. 9.6 of Ref. |119j . 
pp. 182-186 of Ref. [US, pp. 499-501 of Ref. [125], and pp. 180-186 of Ref. [122] . 
Recently, this came again into focus [1711115) . The dynamic algebra so (2d) was 
also discussed as symmetry of Hamiltonians in Rcfs. 11261 1127] . 

The polynomial growth (in d) of the algebra in Theorem [31] to the dynamic 
system of Eqns. ([3"3"]) was identified as the reason why efficient classical simula- 
tion of quadratic, fermionic quantum systems is possible (cp., e.g., pp. 9-10 of 
Ref. [128] . pp. 5-6 of Ref. [HI], and pp. 3 of Ref. [HD]). — Setting n := d- 1 in 
Proposition [27] we obtain 

Corollary 32. Heisenberg-XX chains of d — 1 spin-^ qubits (d > 3 J and two 

individually locally controllable qubits, one at each end, have the dynamic algebra 
so (2d) and can thus simulate a general fermionic quadratic Hamiltonian on d 
levels and vice versa. 

By the second and third statement of Theorem |3"T1 the embedding of so (2d) 
into su(2 d ) splits into two irreducible representations of equal dimension, and 
thus so (2d) acts simultaneously on both components. Therefore, this embedding 
is — besides a doubling of the dimension — equivalent to the embedding of so (2d) 
into su 

r 2 (d-l)i 

via Proposition [27] as readily verifiable by resorting to the Pauli 
basis for so (2d) given there. Referring to Tab. [2] we further remark that the 
embedding of so(2d) into su[2( d_1 )] can arise from a symplectic representation 
(for d = 4fc + 2), an orthogonal one (for d = 4fc), or a unitary one (for d odd). 

Now we illustrate that a controlled spin chain can simulate a quadratic 
fermionic system, while the converse does not hold. To this end, consider the 
case where the general quadratic ('physical') Hamiltonian is supplemented by 
the ('unphysical') linear terms J2p=iJpfp + ^p/p with j p , k p € C. Hermiticity 
implies j p = k p . Assume again that the coefficients j p and k p are real, i.e. j p = j p 
and k p = k p to obtain j p = j p = k p . Thus the linear terms can be written as 

S p =iip(/p + /p)i they are mapped via the Jordan- Wigncr transformation to 
the operators 

d 

H 2 = Y^jp Z^ZXW . (34) 

P =1 p-l d-p 

As will be shown next, this determines the dynamic algebra of a fictitious Hamil- 
tonian system containing quadratic and linear terms (see also the Pauli basis for 
so (2d + 1) given in the proof to Proposition [26]) : 

Theorem 33. Let j p G K (1 < p < d) and the entries of the real antisymmetric 
matrix A and the real symmetric matrix B denote the control functions in the 
Hamiltonian H + Hi, where H and H2 are given by Eqn. (|33p and Eqn. (1341) . 
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respectively. The dynamic algebra so(2d+ 1) of the corresponding control system 
is irreducibly embedded in su(2 d ). 

Proof. Computing commutators of generators — |Z • • • ZXI • • • I from Eqn. 
with generators — |I • • • IZI • • • I from Eqn. (|33a[) , we obtain the additional gener- 
ators — |Z • • • ZYI • • • I. Furthermore, the generators — |I • • • IXZ • • • ZXI • • • I and 
— |I - • • IYZ • • • ZYI • • • I arise from linear combinations of Eqns. P3b[) - (|33cp and 
computing commutators with generators — ^1 • • • IZI • • • I from Eqn. Q33ap re- 
veals the generators — |l - • IXZ • ■ ■ ZYI - ■ ■ I and — |l • • • IYZ ■ ■ • ZXI ■ I. Now 
the Theorem follows by comparing all the generators with the table in the proof 
of Proposition 1261 ■ 

Moreover, setting n := d in Proposition [25] we get 

Corollary 34. Heisenberg-XX chains of d spin-^ qubits (d > 1) and a single 
locally controllable qubit at one end have the dynamic algebra so(2d+ 1) and can 
simulate a general fermionic quadratic Hamiltonian on d levels with its dynamic 
algebra so(2d) ; but obviously not vice versa 5 . 



10.4- Discussion. We analyse three further cases of fermionic Hamiltonians. 
First, consider quadratic Hamiltonians (without linear terms) which are particle- 
number preserving, i.e. A = in Eqn. ([3"3l . Assuming the elements of B are 
control functions, we obtain u(e?) as dynamic algebra (cp. p. 501 of Ref. |125j ). 
Second, the diagonal normal form (see, e.g., App. A of Ref. |112) . Sec. III. 8 of 
Rcf. [Tig, Sec. 3.3 of Ref. PH], and Theorem III of Ref. [EE]) for the Hamil- 
tonian H of Eqn. (gHJ) 

E e p (/t/ p - 1) = E f ~ f ^ (35) 
p—i p—i 

(with E p positive and real) is mapped by the Jordan- Wigncr transformation to 
the G?-qubit operator 

E^l^JZI^I- (36) 

P =1 p-l d-p 

Considering E p (1 < p < d) as controls, we get a (i-dimcnsional abclian Lie 
algebra as dynamic algebra. Third, if we allow for fermionic operators of arbitrary 
order (less than or equal to d), we get su(2 d ) as dynamic algebra. 
In summary, the sequence of dynamic Lie algebras 

su{2 d ) Dso(2d+ 1) D so(2d) D u(d) (37) 

plays a prominent role for fermionic quantum systems as pointed out in Chap- 
ter 22 of Ref. [131] and Sec. IV of Ref. |132) . This sequence of Lie algebras is 
also widely studied in particle physics [13311134] . 



5 as is also illustrated by the unphysical linear terms above 
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10.5. Spinless Hubbard Model with Periodic Boundary Conditions. First, spec- 
ify a spinless version of the Hubbard model (see p. 20 in Ref. |135| or p. 61 
in Ref. 11361 ) with periodic boundary conditions where d + 1 = 1 due to the 
periodicity, d > 2, and t,u EM.: 



H = -tj2(P P f P +i-f,4 + i) 
+«E(4/p-§). 



(38a) 
(38b) 



p=i 



Equation (|38a[) resembles a spinless tight-binding model (see p. 437 in Ref. |135j 
or p. 59 in Ref. [136]) and equals the quadratic Hamiltonian of Eqn. (f29|) with 
A = and 

/0 1 
1 1 
1 • 



B 



V 



1 

1 1 
10/ 



The Jordan- Wigner transform of the Hamiltonian H of Eqn. 



takes the form 



H 



l[{J2b^J YY b^J + i_^_ixx i_^_i ) + (xz_^_zx + yz_^_zy) 

P =1 p-1 d-l-p p-1 d-l-p d-2 d-2 

(39a) 

d 

+ u^I_^_IZI_^_I . (39b) 

P =1 p-1 d-p 

Extending Eqn. (|38bl) to J2 P =i u p(fpfp~ \ ) sucn that it contains site-dependent 
control functions u p E ffi, we obtain (by building on App. B of Ref. [115] ) 



Lemma 35. The dynamic control system corresponding to the Hamiltonian 

d d 



p— 1 p=i 



(40) 



/ias u(d) as dynamic Lie algebra assuming that n p el are controls and d > 2. 

Proof. Let denote the dynamic algebra of the control system. We obtain from 
Eqn. (|39a[) one generator a\ and from Eqn. (|39b[) the generators (0 < p < d) 

Zp ;= ^}^p = ^I***IZl***l. 



p-i 



d—p 
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One can verify on the generators that the d-qubit operators — ^Z---Z and 

— f X)p=i are both elements of the centraliser of td- By comparison with 
Theorem [2U we obtain that td Q so(2e£). As the centraliser in Theorem I3T1 is 
one-dimensional and the centraliser of td is at least two-dimensional, it follows 
that td 5 so (2d). We remark that u(d) is a maximal subalgcbra of so (2d) and 
that su(q) © u(d — q) is a maximal subalgebra of su(d) (see, e.g., p. 219 of 
Ref. |123j ). In particular, u(q) © u(d — q) is a maximal subalgcbra of u(d). The 
Theorem can be directly verified for d € {2,3,4,5}. We assume by induction 
that the Theorem is true for all td' with d! < d. The Theorem follows by in- 
duction if we show that td 2 u (<?) © — <?) holds for any q. We compute the 
commutators a 2 ■= [zi,[z 2 ,ai]] = — |(XiX 2 + YiY 2 ) and a 3 := [z 3 ,[z 2 ,ai]] = 

— |(X 2 X 3 + Y 2 Y 3 ) [using again the notation of Eqn. ([25]!]. For 3 < j < d— 1 we 
have gj := [zj+i, [zj,ai]] = — |(XjX J+ i + YjY J+ i) and by linear combinations 
we obtain the (i-qubit operator := — |(XZ • • • ZX + YZ • • • ZY). We further 
compute the commutators a.5 := [z 2 , [a 2 , 04]] = — ^(IXZZ • • • ZX + IYZZ • • • ZY) 
and a 6 := [z 3 , [03, a 5 ]] = -|(ILXZ • • • ZX + IIYZ • • • ZY). Using the elements gj 
and ae we can build the element a^— §(X)j=3 XjXj+x +YjY J , +i) which together 
with a 2 and the elements Zj generates u(2) © u(d — 2). As a 3 is not contained in 
u(2) © u(d - 2) we proved that t d 2 u(2) © u(d - 2). ■ 



iVofe on i/ie Hubbard Model with Periodic Boundary Conditions and Spin. 
Including the spin a = ± in the Hubbard model gives 



H = — t 



E E(/ 3 

<7 — ± p—1 



t 

p,a 



fp+l,a fp,afp+l cr) 



+ u p(fp, + fp,+ g)(/p,-/pi- 2 )' 



(41a) 



(41b) 



where the anticommutation relations of Eq. ([27)) still hold among operators with 
egtta/ spin values, while operators with different spin values anticommute. The 
spin degrees of freedom just split each of the original levels p into two sub-levels. 
Thus the image of the Hamiltonian form Eq. (|4"Tj) under the Jordan- Wigner 
transformation operates on a space of squared dimension as compared to the 
case without spin and the dynamic algebra is embedded in su(2 2<i ). The drift 
Hamiltonian of Eq. (|41a[) is mapped to Aq © 1 + 1 © Aq with 



d-1 

A) := fVl - IYYI-- 



1=1 p-i 



I 



T_^_IXXI_^_I J + [XZ_-ZX 

p—1 d— 1— p d— 2 c 

The control Hamltonians of Eq. (|41bp are mapped to A p © A p where 



YZ-ZY 



d-l-i 



A,, 



I...IZI-I 

p—1 d—p 



For d = 2, direct computation using the computer algebra system MAGMA 87 
gives the system Lie algebra su(2)©su(2)©u(l) embedded in su(2 4 ). The general 
case appears more intricate and goes beyond the scope of this work. 
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11. Bosonic Quantum Systems 

Finally we comment on the bosonic case. As opposed to the Pauli principle in 
the fermionic case, in bosons the occupation number n p is not bounded and — 
even for a finite number d of levels — the dynamic algebra of Hamiltonians of 
arbitrary order need not be finite unless the particle number is also bounded. 
Yet, the dynamic algebra for quadratic (pair-interaction) Hamiltonians is given 
by the real symplectic algebra sp(2d, K) (see, e.g., p. 36 of Ref. |114j . p. 186 of 
Ref. |120j . or p. 501 of Ref. |125j ). We have not yet found an appropriate spin 
system that would be dynamically equivalent to the compact real form sp(d) of 
a quadratic bosonic system with algebra sp(2d, K). However, in Secs. l9.2H9.3l we 
have already presented spin systems with dynamic algebras sp(2" _1 ) which are 
actually more powerful than required and contain the compact real form sp(d) 
of a quadratic bosonic system with algebra sp(2d, R). For further analysis of 
the bosonic case, the Holstein-Primakoff transformation may be of help (see, 
e.g., p. 78 of Ref. |136j ). — Finally, the results of mutually simulating quantum 
systems are summarised in Tab. [3] 



12. Outlook: 

Quantum Simulation as an Observed Optimal-Control Problem 

Clearly, in view of experimental settings, one may take a more specific point 
of view by comparing the time course of two observed bilinear control systems 
(Sfi), p — a,b with respect to (i) a set of Hcrmitian (and mutually orthogonal) 

observables and C^) with v, v' elC {1,2,..., iV 2 — 1}, (ii) the initial states 
Pq, (iii) a given time interval [0, T], and (iv) admissible controls Uj(t) € C R 

p»(t) = -i + J2 > nt)] with p"(0) = p% (42) 

(C)Z(t) = tr{(C£) V(i)} with {Ca .61. (43) 

Now the comparison resorts to the expectation values (C)^(t) via states p^(t) : 
drifts Hq , controls H^, and control amplitudes u^(t). Note that {C^'} and 
{^V } ncc d not coincide, but if S a shall simulate Si, it is convenient to re- 
quire |{C^}| > |{C^}| so that (by invoking the above orthogonality of the 
observables with respect to the Hilbert-Schmidt scalar product) one can ensure: 

rank spanR-(C,^} > rank spani{{C^}. 

Now for simultaneous measurement, it is useful to pick several observables 
C£f as long as they are compatible (mutually commute), or, more generally, as 
long as they are mutually non-disturbing in the sense of the recent findings 
in Ref. [23j . Simultaneous expectation values are conveniently collected in the 
observation vectors 

[(CrW^KC)^),^)^),...}* with p = a,b. (44) 

Likewise, we define the respective dynamic system algebras of S a and Sb as 

e M := (iHg,iH?\j = l,2,...,m^) Lio with y, = a, b . (45) 
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Table 3. Summarising Overview on Simulating Quantum Systems 
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n by the Jordan- Wignor 



Clearly, t a D £ b implies S a can simulate However, if S a comes with a larger 

set of observables {C^}, the above condition is still sufficient, but it is no longer 
necessary. This is analogous to the fact that in quantum systems controllability 
implies observability, whereas the converse need not hold [35] (for details see 
|37j). In classical systems, however, controllability and observability are dual to 
one another (see, e.g., [137) ). since no observables accounting for the quantum- 
specific measurement process are involved. — Now the notion of weak simulation, 
for which simulability can be seen as a strong condition, comes naturally: 

Proposition 36. A dynamic system S a can weakly simulate another dynamic 
system Sb in time interval [0,T] and with respect to the two sets of observables 
{C"} and {C^}, if there exists a pair of initial conditions and p$ (reachable 
form the respective equilibrium states) and two sets of admissible control vectors 
u a {t) and u^t') such that M[(C) a (t)} = [(C) b (t')} for all t G [0,T] and t' £ 
[t(0),t(T)], where r(t) is a bijective function oft for all t £ [0,T] and M is a 
map M : R n — > R"\ [(C) a (t.)} h-> [(C) b (t)} with n>m. 

As will be described elsewhere, the previous proposition motivates to view 
simulability as a generic precondition to formulate weak quantum simulation 
as an optimal-control task: minimise | \AI [(C) a (t)} — [(C) 6 (<')]||| subject to the 
differential equations of motion given in Eqn. (|42l) . 



13. Conclusion 

Often the presence or absence of symmetries in quantum hardware architectures 
can already be assessed by inspection. Given the system Hamiltonian as well as 
the control Hamiltonians, (i) we have provided a single necessary and sufficient 
symmetry condition ensuring full controllability, and (ii) in view of practical 
applications we have shown easy means (solving systems of homogeneous linear 
equations) to determine the symmetry of the dynamic system algebra { merely in 
terms of its commutant or centraliser 6'. If the system Hamiltonian corresponds 
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to a connected coupling graph, the absence of any symmetry can be further 
exploited to decide full controllability: it means the dynamic system algebra 
is irreducible and simple. Now conjugation to simple orthogonal or symplectic 
candidate subalgebras can again be decided solely on the basis of solving systems 
of homogeneous linear equations. The final identification task can now be settled 
because here we have given a complete list of irreducible simple subalgebras 
of su(.ZV) compatible with the physical constituents as a dynamic pseudo-spin 
system. This avoids the usual and significantly more costly way of explicitly 
calculating Lie closures. We have thus made precise and easily accessible the 
following four conditions ensuring full controllability of a dynamic qubit system 
in terms of its system algebra { C su(N) : 

(1) the system must not show a symmetry (t must have a trivial centraliser £'), 

(2) the coupling graph of the control system must be connected, 

(3) the system algebra t must not be given in a symplectic or an orthogonal 
representation, and finally 

(4) if 6 is given in a unitary representation, it must not be on the list of 
proper irreducible unitary simple subalgebras of su(N), in particular, 6 ^ te- 

The system algebra completely determines the possible dynamics of controlled 
Hamiltonian systems. Therefore, the lattice of irreducible simple subalgebras to 
su(N) given here also provides an easy means to assess not only the somewhat 
easier cases of mutual simulability but also the more intricate cases of simulability 
with least overhead of dynamic systems of spin or fermionic or bosonic nature. 
In a number of examples (see also Tab. [3]), we have illustrated how controlled 
quadratic fermion and boson systems can be simulated by spin chains and in 
certain cases also vice versa. 

Finally, since full controllability entails observability (while in the quantum 
domain the converse does not necessarily hold), symmetry constraints immedi- 
ately pertain to observability as discussed in more detail in Ref. (37j . 
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APPENDIX 

A. Tensor-Product Structure in Qudit Systems with Many-Body 
Interactions 

For quantum simulation, we generalise the discussion such as to embrace qudit 
systems with (effective) many-body interactions. Treating them as control sys- 
tems embedded in su(N), now su(di) ffisu^) ffi • • • ®su(d n ) is a tensor-product 
structure of su(N), where Yl" = idj = N and dj > 2. We consider the sub- 
algebras su(dj) as subsystems of the tensor-product structure. We say that the 
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tensor-product structure f)i ffi f)2 ffi • • ■ © fin is a refinement of the tensor-product 
structure su(di) ©5u(d2) © • • • ©5u(c?„) if f)j is either equal to su(dj) or equal 
to su(cj ; i) ffisu(cj j2 ) © ••• ®su(cj !mj ), where Ofc^i c j',fe = dj and rrij > 2. We 
call ()i © f)2 © • • • © f) n a proper refinement if there is one j such that f)j 7^ 
su(dj). For a given quantum system in su(N), there exists a common refinement 
su(pi) ffisu(p2) © • • • ffisu(p„) of all tensor-product structures, where Ylj = i Pj is 
the factorization of N into prime numbers. The common refinement is unique 
up to permutations of subsystems. 

Now with respect to tensor-product structure su(di) ©su(d 2 ) © • • • ffisu(c?„), 
again we write Hamiltonians as a linear combination (cfc G R) 

m 

H = J2 ckn k (46) 
fc=i 

of elements Hk = — fCHfc.i © %fc,2 © • • -®Hk,n) forming a tensor basis of su(N). 
The elements S 2% U are choosen relative to bases 

{-tA\Ae B e :={B /i i J B /i2 ,...,B <i((1<) 2_ 1 }} 

of su(de). In addition, we assume that the order 

ord(?4) :=#{£: H k ,t + W > 1. 

Recall that the Hamiltonian H has a coupling graph, if its order ovd(H) = 
ord(^feLi CkHk) '■= max({ord('Hfe) \k = 1, . . . , m}) is equal to 2, which is the 
case in pairwise coupling interactions. The vertices j are given by the subsystems 
su(dj) and we get an edge between the nodes k\ and fc 2 with k\ 7^ fc 2 if there 
exists a in Eqn. (T4"6")) such that {fci,fc 2 } = {j : 'Hfcj ^ 1^-}- If all control 
Hamiltonians are local, i.e. are contained in su(di) ©su(c?2) © • • • ©su(g?„), then 
we say that the coupling graph of the drift Hamiltonian Hd is the coupling graph 
of the control system. 

We say a control system on su(N) is weakly connected, if the dynamic algebra 
£ contains for each proper partition of its tensor- product structure in (X\ U I2 — 
{l,2,...,m},I 1 nl 2 = {}) 

f)i = ® j e xiSU(dj) and f) 2 = © j e x 2 su(dj) 

an element of su(N) \ [t)iffif) 2 ]. For Hamiltonians H of ord(-ff) = 2, this is 
equivalent to the fact that the coupling graph is connected. We will also use the 
stronger notion of a directly connected control system for which the dynamic 
algebra t contains an element of su(dj 1 dj 2 ) \ [su(d :)1 ) ®su(dj 2 )] for each pair of 
subsystems su(d J1 ) and su(dj 2 ) with j\ 7^ j2- — With these notions, Theorem [1] 
generalises as follows. 

Theorem 37. Consider a bilinear control system onsu(J}" =1 dj), where dj > 2. 
Assume that the subsystems Six(dj) with j G {1, . . . , n} are independently fully 
controllable so the dynamic algebra t 3 su(di) ©su(d 2 ) © • • • ®su(d n ). The con- 
trol system is fully controllable, i.e. t = su(J}" =1 dj), if and only if the control 
system is directly connected. In particular, { = Su(JX, = i dj) is simple. 
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Proof. The 'only if'-direction is obvious. We prove the 'if -direction. First, we 
assume that n = 2. As the subsystems are independently fully controllable, 
we obtain t D su(di) © su(d,2). The dynamic algebra 6 contains an element of 
su(did,2)\[su(di) ©su(d 2 )], as the control system is directly connected. It follows 
from Thm. 1.3 of [55] that su(di) ©511(^2) is a maximal subalgebra of 5u.(d\d2)- 
As t 2 stl(di) ffisu(d 2 ), this proves { = 511(^1^2)- The general case follows by 
induction on the number of subsystems. We remark that su(]^["_ 1 dj) is simple 
so the last assertion follows. ■ 

This complements results on the controllability of quantum circuits |138j . 
where the controllability of continuous and discrete sets of unitary transforma- 
tions is considered. In particular, Theorems 4.1 and 4.2 of Ref. |139j (see also 
140 ) rely also on the maximality of the subgroup of local operations on two 
qudits [i.e. on SU(<i 2 ) 2> SU(d) © SU(d)]. Our controllability proof can be com- 
pared to proofs relying on Cartan decompositions (see Thm. 5 of Ref. [141| and 
Prop. 2.4 of Ref. [H]). Unfortunately, one cannot substitute 'directly connected' 
with 'weakly connected' in Theorem 1371 

Example 38. Consider a bilinear control system on su(8) with the tensor-product 
structure su(2) ffisu(2) ffisu(2). We assume that the subsystems are indepen- 
dently fully controllable, i.e. 

D (iXII, iYII, iZII, ilXI, TYI, TZI, TLX, TIY, TIZ) Lic . 

In addition, we have a drift Hamiltonian Hd = ZZZ. The control system is weakly 
connected and { acts irreducibly. The dynamic algebra is t = sp(4) ^ su(8) and 
hence the system is not fully controllable. 

B. Connected Control Systems in Qudit Systems with Many-Body 
Interactions 

In this Appendix we build on Sec. [S] and discuss a more general notion of con- 
nected control systems in qudit systems with many-body interactions which do 
not necessarily have a natural coupling graph. We freely use the notation of 
Appendix [Al 

Recall Example [S] of Sec. [5] Motivated by this example one might conjecture 
that the dynamic algebra I is simple if the control system is weakly connected 
and t acts irreducibly. Unfortunately, this is not true. 

Example 39. Assume that we have a bilinear control system on su(8) with two 
subsystems corresponding to the tensor-product structure su(4) © su(2) . On the 
first subsystem we pick f)i = (iXII, iYII, iZII, ilXI, ilYI, TZI)Li C as the local dy- 
namic Lie algebra. On the second subsystem we pick the local dynamic Lie alge- 
bra f)2 = (TIX, TIY, TIZ)Lio- In addition, wc have a drift Hamiltonian Hd = IZZ. 
The control system is weakly connected and t acts irreducibly. We obtain that 
the dynamic Lie algebra is t = su(2) ©su(4). It is neither simple and nor fully 
controllable. In particular, the dynamic Lie algebra does not respect our chosen 
tensor-product structure. The problem is that the control system 
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is not weakly connected w.r.t. the tensor-product structures su(2)ffisu(4) and 
su(2)©su(2)§su(2). 

Generalising Sec. [51 we say that a control system is connected, if the dynamic 
algebra 6 contains an element of su(N) \ [su(ei) ©su(e2)] for each tensor-product 
structure su(e±) ©su(e2) with TV = eie2 and ex, C2 > 2. For control systems with 
pair interactions this definition is equivalent to the one given in Sec. [5J 

Lemma 40. The following are equivalent: 

(1) The control system is connected. 

(2) The control system is weakly connected w.r.t. the common unique refinement 
of its tensor-product structure. 

(3) The control system is weakly connected w.r.t. any tensor-product structure. 

■ 

We now generalise Theorem [5] and prove that the dynamic algebra £ is simple 
if its centraliser is trivial and the corresponding control system is connected. 

Theorem 41. Assume that the dynamic algebra t of a bilinear control system 
on su(N) has a trivial centraliser ¥ . Then one finds: 

(1) The dynamic algebra 6 is given in an irreducible representation. 

(2) If I is semi-simple but not simple, then t ^ su(N) and the control system is 
not fully controllable. 

(3) The dynamic algebra 6 is simple iff the control system is connected. 

Proof. (1) immediately follows from t' being trivial and Lemma [31 while (2) is 
obvious, as su(N) is simple. We now prove the 'if '-part of (3). We obtain from 
Lemma 0] that t is simple or semi-simple. In the following, we assume that 4 
is not simple. Thus 6 is a irreducible semi-simple (but not simple) subalgebra 
of su{N). Using Thm. 2.1 of Ref. [88], it follows that 4 = ti§t 2 § ••• ©« m 
and that the 6, are irreducible simple subalgebras of some su(dj) such that 
t C su(di) ffisu(d2) © • • • ffisu(d m ), YTjLi dj = ^> an d m > 2. In particular, we 
can choose two non-zero algebras 

hi = ffijeiifij and fo 2 = ©jex 2 %, 

where t = fhffif) 2 C su(ci) ffisu(c 2 ), X x U X 2 = {l,2,...,m}, X 1 n T 2 = {}, 
and c\C2 = N. As the control system is connected, the dynamic algebra 6 con- 
tains an element of su(N) \ [su(ci) ©su(c2)] for each tensor-product structure 
su(ci) ©su(c2). This is a contradiction to t C su(ci) © su(c2) and the 'if'-part 
of (3) follows. To prove the 'only if'-part of (3) we assume that the control sys- 
tem is not connected. It immediately follows that the dynamic algebra has to 
be a (non-trivial) direct sum. Thus it cannot be simple, which proves the 'only 
if'-part by contradiction. ■ 

In important special cases more convenient conditions hold: 

Corollary 42. Given a bilinear control system on su(N), where the centraliser 
t' of the dynamic algebra t is trivial. We obtain: 

(1) Assume that the subsystems of the tensor-product structure are independently 
fully controllable. The dynamic algebra t is simple if and only if the control 
system is weakly connected. 
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(2) Assume that the tensor-product structure of the control system is given by 
su(pi) 0su(p2) © • • ■ ®su(p n ), where YYj=iPj * s a factorization of N into prime 
numbers. For example, pj = 2 for all j . The following are equivalent: 

(a) The dynamic algebra £ is simple. 

(b) The control system is weakly connected. 

(c) The control system is connected. 

Proof. We first prove the case of (1). As the subsystems are independently 
fully controllable, any irreducible semi-simple (but not simple) dynamic algebra 
i D su(di) ®su(d2^ © •• • ©su(<T m ) has to be (irreducibly) contained in the al- 
gebra f) = su(d[) ®su(d' 2 ) © • • • ®su(d' m ) where su(di) ffisu(d 2 ) © • • • §su(d m ) 
is a refinement of the tensor-product structure f). All these cases are excluded 
as the control system is weakly connected, and (1) follows along the same lines 
as Theorem HU As su(pi) ©su(p2) © • • • ®su(p n ) is the common unique refine- 
ment of all tensor-product structures in the case of (2), the control system is 
weakly connected if and only if it is connected (by Lemma |30j) and (2) follows 
by Theorem EH ■ 

C. Computational Techniques for Representation Theory 

For computationally exploiting Lie theory to list all irreducible representations 
of a given dimension TV for all irreducible simple subalgebras of su(TV), a self- 
consistent frame is indispensable. It requires the highest weights and the di- 
mensions of their respective representations to be linked to the classification by 
the standard Dynkin diagrams. Here we explicitly give all the details in such 
a consistent frame, since combining different literature sources runs the risk of 
arriving at erroneous results due to possibly inconsistent conventions. 

In particular, the appendix is meant to complement Sec. [5] It describes the 
methods we used to compute the irreducible simple subalgebras of su(TV) and 
their inclusion relations. 



C.l. Highest Weights and Dimension Formulas. The irreducible simple subal- 
gebras of su(TV) are found by enumerating for all simple Lie algebras all their 
irreducible representations of dimension TV. The irreducible representations can 
be enumerated using highest weights {x\, . . . ,xe) which are (non- negative) inte- 
ger vectors. The length I of the highest weight is given by the rank (i.e. dimension 
of the maximal abelian subalgcbra) of the considered Lie algebra. Details on the 
theory of highest weights can be found in Chap. IX, Sec. 7 of Ref. [76] . 

Different orderings for the coefficients Xi of the highest weights are used in the 
literature. We use the so-called Bourbaki ordering which is detailed in Tab.@]by 
numbering the nodes of the Dynkin diagrams (see Chap. VI, Sec. 4.2, Thm. 3 of 
Ref. [75] ) for the compact simple Lie algebras. In Tab. [5] we present the highest 
weights and dimensions for the standard representation of each compact simple 
Lie algebra. We put highest weights together if they differ only w.r.t. an outer 
automorphism, i.e. an permutation which leaves the Dynkin diagram invariant. 
Note that the standard representation is the lowest-dimensional (non-trivial) 
representation [with the exception of so(3), so(5), and so(6)] and is typically 
used to introduce the corresponding Lie algebra in matrix form. 
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Table 4. The Compact Simple Lie Algebras and their Dynkin Diagrams 
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Wc already remarked in Sec. [6] that the dimensions of irreducible representa- 
tions can be efficiently computed using computer algebra systems such as Li E |86j 
and MAGMA [57] via Weyl's dimension formula. Now we present explicit formu- 
las for the dimensions, which allowed us to speed up the computation of the 
dimensions considerably. While for su(£ + 1), so(2£ + 1), and sp(£) these formu- 
las can readily be found on pp. 340-341 of Ref. [72], we had to correct the one 
for so(2£), since we could not find a reference with the proper formula either. 

Lemma 43 (Classical Lie algebras). Given the highest weight (xi, . . . , xg) 
the dimensions of the corresponding irreducible representations are: 



1. su(£+l) : dim = J]i 



<i<j<£+ 



1 I i~* i 
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Table 5. The Compact Simple Lie Algebras and their Standard Representations (with the 
corresponding dimensions) 



algebra highest weight(s) dim 



(0,...,0,1) 



so(3) 


(2) 




3 


so(2£+l) 
l>2 


ObJV 
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..,0) 
-1 
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sp(£) 


(1,0^ 

i 


..,o) 

-1 


21 


so(2) a 


(1) 




2 


so(4,) a 


(1,1) 




4 



"not simple 



algebra highest weight (s) dim 



so (6) (1,0,0) 



so(8) 


(1,0,0,0), 
(0,0,1,0), 
(0,0,0, 1) 


8 


60 (2f) 


(1,0,..., 0) 
e-i 


21 


e6 


(1,0,0,0,0,0), 
(0,0,0,0,0,1) 


27 




(0,0,0,0,0,0,1) 


56 




(0,0,0,0,0,0,0,1) 


248 


f4 


(0,0,0, 1) 


26 


02 


(1,0) 


7 



2. 5 0{2l+ 1) : dim = JJ, , ,. . {l + Si±^±^fe±^^ | 
x II: , . {l + x II, , , {l + 2(x ^:^ +xe } 

3. SP (£) : dim = H, , ,. , {l + Xi+ ~ j +*>- 1 } X H, , , {l + 

tt f-, | Zi + ---+a;j-i+2(a:j + ---+3;f) ^ 

A lll<i<j<£ \ ' 2<+2-t-j J 

I *or>r : : dim = H : , , , {l + } X H, , . , {l + Ei±^±2i | 

T-r fi , hrj_i+2(xjH hri_ 2 )+3:i-i+:ci 1 ■ 

X lll<i<j<^-l I 1 "I 2£-i-j J 

Here we present the dimension formulas for the exceptional Lie algebras only 
for 02 (cp. Ref. |142j . pp. 257-258) and f4, ommitting the even longer and more 
complicated ones for Zq, e^, and es. We remark that these formulas are — in 
principle — well known but are usually not given in the literature due to their 
complexity. 

Lemma 44 (92 and j^). Given the highest weight (xi, . . . , xi) the dimensions 
of the corresponding irreducible representations are: 



2:ri+3x2 ' 
5 



1. 02 : dim = (1 + a? a )(l + a?i) (1+^4^) (1+^^) (l+^f^) (l + 

2. f 4 : dim = (1 + x 4 )(l + x 3 )(l + z a )(l + a?i) (l+ 2a ^ 4 ) (l+^^) 

x ^1 | a 2+a3+XA ^ ^ i 2x 2 +X 3 j ^ [ Xi+X 2 +X 3 j ^ j 2x 2 +X 3 +X 4 ^ ^ 1 X1+X2+I3+I4 ^ 

x ^ 1 Xi+2x 2 +X 3 \ /-^ 1 2x 2 +2x 3 +X 4 ~j ^ 1 Xi+2x 2 +X 3 +X 4 \ /-^ 1 2xi+2x 2 +x 3 ^ 

x ^ [ 2'i+2x 2 +2x 3 +X4 ^ A | 2xi+2x 2 +x 3 +x 4 ^ ^ [ xi+3x 2 +2x 3 +x 4 S 

x ^ [ 2xi+2x 2 +2x 3 +x 4 ^ j 2xi+3x 2 +2x 3 +x 4 \ ^ | 2xi+4x 2 +2x 3 +x 4 ^ 

x ^ j 2xi+4x 2 +3x 3 +x 4 j ^ [ 2xi+4x 2 +3x 3 +2x 4 ^ 
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Proof. Computational explicit dimension formulas for the exceptional Lie alge- 
bras were obtained using the computer algebra system MAGMA [87] via Weyl's 
dimension formula. I 

We emphasize that in order to compute the dimensions efficiently, one has 
to use the dimension formulas in the given factorized form. That is, one has to 
evaluate each factor and multiply the results. The alternative of evaluating the 
multiplied formula is considerably less efficient. 



C.2. Enumerating Representations. The aim of determining the irreducible sim- 
ple subalgebras of su(N) for a given N is reached by enumerating for all simple 
Lie algebras all their irreducible representations of dimension N. Therefore, we 
have to enumerate for all simple Lie algebras all highest weights (xi, . . . , xi) 
corresponding to irreducible representations of dimension N . In doing so, how 
can one reduce the combined search space of Lie algebras and highest weights? 

To this end, recall that the standard representation is the lowest-dimensional 
(non-trivial) representation [with the exception of so(3), so(5), and so(6)]. It 
follows from the dimension formulas for the standard representations in Tab. [5] 
that only a finite number of Lie algebras have irreducible representations of 
dimension equal (or less than or equal) to a given N . Thus we have to search 
only through a finite set of Lie algebras. In addition, we have to consider merely 
one instance of isomorphic Lie algebras [su(2) = so(3) = sp(l), so(5) = sp(2), 
and su(4) = so (6)] and can neglect so (2) and so (4) as they are not simple. It 
follows from Chap. IX, Sec. 8.5, Cor. 2 of Ref. [73] that for a Lie algebras the 
set of irreducible representations of dimension less than or equal to N is finite. 
We obtain: 

Lemma 45. Each Lie algebra can only have a finite number of irreducible repre- 
sentations of dimension equal (or less than or equal) to a given N . Furthermore, 
only a finite number of Lie algebras have any irreducible representations of di- 
mension equal (or less than or equal) to a given N . M 

It follows from Lemma US] that the search space for the highest weights is finite. 
Using the following Lemma, one can obtain stopping criteria for the search for 
the highest weights with dimension equal (or less than or equal) to a given N. 

Lemma 46. For a given Lie algebra, let dim [a;] denote the dimension of an 
irreducible representation with highest weight x = (x\, . . . , xi). 

1. The dimension is strongly monotonic ascending in each entry Xi of the highest 
weight: dim[(xi, . . . , Xi + 1, . . . , xi)] > dim[(xi, . . . , x,, . . . , xi)]. 

2. Let e 1 denote the vector such that (e*)j = 8ij. If V, Xj > 1, then 

dim[x] > min[{max[dim(e t ), dim(e- 7 )]}i^j U {dim(2e l )}i<,<^]. 

Proof. Sec, e.g., Cor. 5.2 and Cor. 5.4 of Ref. [H3]. ■ 

Let us fix a Lie algebra. We start our search for the highest weights of dimen- 
sion less than or equal to N with all vectors x = (xi, . . . , xi) such that x i = 1- 
In each step, we compute the dimension corresponding to the highest weight 
x = (x\ , . . . , xi) . If dim[x] < N, we include x into the list of highest weights with 
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dimension dim Ixl and wc branch our search to all + l,...,x t ). If 

dim[x] > N, we prune this branch in our search tree (Part 1 of Lemma |46| . We 
have to search through all Lie algebras such that the lowest-dimensional (non- 
trivial) irreducible representation is less than or equal to N. One can further 
reduce the search space with respect to potential Lie algebras by using knowl- 
edge on the sccond-lowcst-dimcnsional (non-trivial) irreducible representations: 

Theorem 47. For a given Lie algebra, let y = . . . , yi) denote the highest 
weight of the second-lowest-dimensional (non-trivial) irreducible representation. 

1. For su{£+l) and £ > 3, we obtain y = (0, 1, 0, . . . , 0) or y = (0, . . . ,0,1,0). 
In addition, dim[y] = £(£ + l)/2. 

2. For so(2£+l) and £ > 7, we obtain y = (0, 1, 0, . . . , 0) and dim [y] = (2£+l)£. 

3. For sp(£) and I > 4, we obtain y = (0, 1, 0, . . . , 0) and dim[y] = 2£ 2 - £ - 1. 

4. For so{2£) and £ > 8, we obtain y = (0, 1, 0, . . . , 0) and dim[y] = 2£ 2 - I. 

Proof. Again, let e z denote the vector such that {e l )j = Sij. We apply Lemma l46l 
in each of the following cases: 

1. su(£+l): Recall that dim[e r ] = ( t + 1 ) for 1 < r < £. It follows that dim[e 2 ] = 
£(£ + l)/2 for £ > 3. One can deduce from Lemma |4"31 that dim[2e 1 ] = (£ + 
l)(£ + 2)/2 and that dim[(l, 0, . . . , 0, 1)] = £{£ + 2) for I > 2. We obtain that 
dim[e 2 ] < dimpe 1 ] < dim[(l, 0, . . . , 0, 1)] for £ > 3. The first part follows. 

2. so(2£ + 1): Recall that dim[e f ] = 2 £ and dim[e r ] = ( 2e + 1 ) for 1 < r < £ - 1 
(see, e.g., p. 340 of Ref. [72]). lt follows that dim[e 2 ] = (2£ + l)£ for £ > 3. 
We obtain that dim[e 2 ] < dim[e f ] for £ > 7. One can deduce from Lemma 1431 
that dim[2e 1 ] = £(2£ + 3) > dim[e 2 ] for t > 3. The second part follows. 

3. sp{£): Recall that dimfe 1 ] = 2£ and dim[e r ] = ( 2 /) - (^J for 2 < r < £ (see, 
e.g., p. 341 of Ref. [72]). It follows that dim[e 2 ] =2£ 2 -£-l. Wc obtain that 
dim[e r ] = 2 %% 2r ( 2f + 2 ). One can deduce that dim[e 3 ] = |£ 3 - 2£ 2 - \£ > 
dim[e 2 ] for £ > 4. If r > 4, it follows that dim[e r ] - dim[e 2 ] = 2 %% 2r ( 2£ + 2 ) - 
2£ 2 +£+l > T L [ ( 2e + 2 )-2£ 2 +£ + l = \P - 2£ 2 + \£ + 1 > for £ > 4. One 
can obtain from Lemma H31 that dim[2e 1 ] = 2^ 2 + £ > dim[e 2 ]. The third part 
follows. 

4. so(2£): Recall that dim^" 1 ] = dim[e f ] = 2 f ~ x and dim[e r ] = ( 2 /) for 1 < r < 
£ - 2 (sec, e.g., p. 341 of Ref. [72]). It follows that dim[e 2 ] = 2£ 2 - £ for £ > 4. 
We obtain that dim[e 2 ] < dim[e £ ] for £ > 8. One can deduce from Lemma l43l 
that dim[2e 1 ] = 2£ 2 + £ - I > dim[e 2 ] for £ > 4. The fourth part follows. ■ 

Now one obtains bounds on £ such that the dimension of the second-lowest- 
dimcnsional (non-trivial) irreducible representation is greater than N: 

Corollary 48. For a given Lie algebra, let y — . . . ,yg) denote the highest 
weight of the second-lowest- dimensional (non-trivial) irreducible representation. 

1. For su(£ + 1) and £ > 3, we obtain: dim[y] > N <^ £ > y^TJI+2N - 1/2 

2. For so{2£ + 1) and £ > 7, we obtain: dim[j/] > N £ > (x/TTsT/V - l)/4 

3. For sp{£) and £ > A, we obtain: dim[j/] > N ^ £> (y/9 + 8N + l)/4 

4. For so{2£) and £ > 8, we obtain: dim[y] > N & £ > {y/l + 8N + l)/4 ■ 
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Now wc explain how to use Corollary [?5] in order to reduce the search space. 
Consider su(k + 1) and k > 3. If N > k + 1 but k > ^/T/I+2N - 1/2 then 
the standard representation of su(k + 1) occurs with dimension less than or 
equal to N. But no other (non-trivial) irreducible representation of su(k + 1) 
has dimension less than or equal to N. We include the highest weight of the 
standard representation in the list corresponding to the dimension fc + 1. But we 
do not have to search for other irreducible representations. The search space is 
reduced from a size linear in N to a square-root in N. 



C.3. Inclusion Relations. Once having obtained all irreducible simple subalge- 
bras of su(N), one can determine their inclusion relations following the work of 
Dynkin [55] (see, e.g., Chap. 6, Sec. 3.2 of Ref. [U]). Refer to PH] for related 
literature. For example, Refs. |145U146j generalise the work of Dynkin [55] to 
classical and exceptional Lie algebras over prime fields. References [79111471 .148 
contain most recent findings. It follows from Theorem 1.5 in |88j that almost 
all irreducible simple subalgebras (of dimension dim) are maximal in su(dim), 
sp(dim/2), or so(dim). Relying on Tab. I of [55], the exceptions are listed in 
Tab. [6] which contains irreducible simple subalgebras of su(dim) including the 
algebra in which the subalgebra is maximal. In addition, the highest weights 
of the corresponding representations as well as the type of the subalgebra (s 
for symplectic, o for orthogonal, and u for unitary) are given. For reference, 
we give the Malcev classification JS3] (see also, e.g., 88. 102.55]) of symplectic, 
orthogonal, and unitary representations: 

Theorem 49 (Malcev). Let x = [xx, . . . ,xg) denote the highest weight corre- 
sponding to an irreducible representation <j> of a Lie algebra t. As 4> is irreducible, 
the different cases of symplectic, orthogonal, and unitary representations are mu- 
tually exclusive: 

1. t = 5u(£ + l); 6 

(a) <p is symplectic if x is symmetric, {£ mod 4) = 1, and xiu_ i) div 2)+i is °dd. 

(b) cf> is orthogonal if x is symmetric as well as either (i) (£ mod 4) = 1 and 

div2)+i is even or (ii) (t mod 4) ^ 1. 

(c) <p is unitary if x is not symmetric. 

2. t = so(2£+ 1) for£> 2: 

(a) <p is symplectic if {I mod 4) € {1, 2} and xg is odd. 

(b) </> is orthogonal if either (I mod 4) € {0, 3} or xi is even. 

3. 4 = sp(£) for £ >2: (j) is symplectic if J2i<2j+i<e x 2j+i is odd (j'eNU {0} ). 
Otherwise, </> is orthogonal. 

4. t = so{2£) for£>3: 

(a) <p is symplectic if {£ mod 4) = 2 and xg~i + xe is odd. 

(b) 4> is orthogonal if either (i) {£ mod 4) = 2 and xt-\ + xg is even, 
(ii) {£ mod 4) = 0, or (in) £ is odd and X£-i = xp. 

(c) 4> is unitary if £ is odd and xt-\ ^ xg. 

5. t = Q2, i = f4, or t = eg : <j) is always orthogonal. 

6. £ = t§: (f> is orthogonal if X\ = Xq and X3 = x§. Otherwise, <j) is unitary. 

7. t = e*?: (f> is symplectic if xi + x§ + 27 is odd. Otherwise, </> is orthogonal. M 



6 Recall that div denotes integer division, e.g., (5 div 2) = 2. 
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Table 6. Irreducible Simple Subalgebras not Maximal in su(dim), sp(dim /2), or so(dim) 



subalgebra 


type highest weight (s) 


algebra 


highest weight(s) 


dim 


su(f+l) 

l>4 


u 


(1,0,1,0,... ,0), 
(0,... ,0,1,0,1) su 


W+l)/2] 


(0,1,0,.. .,0), 
(0,... ,0,1,0) 


3(T) 


su(£+l) 

e>3 


u 


(2,i,o„ ..,o), r 

(0,... ,0,1,2) 




(0,1,0,.. .,0), 
(0,.. .,0,1,0) 


H l i 3 ) 


su(2) 


o 


(6) 


02 


(1,0) 


7 


su(6) 


o 


(0,1,0,1,0) 


sp(10) 


(0, 1,0,..., 0) 


180 


so(4/c+3) 

fe>l,m>l 


s/o b 


(0, . . . ,0,m) 


so(4fc+4) 


(0,...,0,m,0), 
(0,...,0,0,m) 


c 


(but not k = m=l) a 










so(9) 


o 


(1,0,0, 1) 


so(16) 


(0,.. .,0,1,0), 
(0,.. .,0,0,1) 


128 




o 


(0,2,0) 


s r( 1 1 


(ni n n n n 

^u, x, u, u, u, u, 


90 


nnf"^ 
sp(3) 


s 


(0, 2, 1) 


S P(< I 


^U, U, 1, u, u, u, 




so(10) 


u 


(0,1,0,1,0) 
(0,1,0,0,1) 


su(16) 


(0,0,1,0,.. .,0) 
(0,... ,0,1,0,0) 


560 


so(12) 


o 


f f\ A A 1 n a\ 

(0, 0, 0, 1, 0, 0) 


sp(16) 


(0, 1,0,0, ...,0) 


495 


so(12) 


s 


(0,0,1,0,1,0) 
(0,0,1,0,0,1) 


sp(16) 


(0,0, 1,0, ... ,0) 


4928 


t6 


u 


(0,0,1,0,0,0) 
(0,0,0,0,1,0) 


su(27) 


(0, 1,0, 0,0,..., 0) 


351 




u 


(0,1,1,0,0,0) 
(0,1,0,0,1,0) 


su(27) 


(0,0,0, 1,0,..., 0) 


17550 


17 


o 


(0,0,0,0,0,1,0) 


5p(28) 


(0, 1,0, ...,0) 


1539 


t7 


s 


(0,0,0,0,1,0,0) 


5p(28) 


(0,0, 1,0, ... ,0) 


27664 


t7 


o 


(0,0,0,1,0,0,0) 


5p(28) 


(0,0,0, 1,0,..., 0) 


365750 


t7 


s 


(0,1,1,0,0,0,0) 


5p(28) 


(0,0, 0,0, 1,0,..., 0) 


3792096 


82 

m>2 


o 


(m,0) 


so(7) 


(m, 0, 0) 2m ~ 





"corrected, for k = m = l wc have s O (7) Cso (8) CSU (8) 
& if (fc + is odd then s else o 

^corrected, { ("*+»• ~ ') ' ("* + m " ') ] 



Results for dimension dim < 16 can be found in Tab. where the irreducible 
simple subalgebras of su(dim) are given again with their type (s for symplectic, 
o for orthogonal, and u for unitary) plus the highest weight of the corresponding 
irreducible representation. This information is essential for deriving Tab. [1] 



C4- Examples. Wc illustrate our methods by two concrete examples: 

Example 50. We use the methods of Appendix IC.2I in the case of dimension 
N = 7 and compute the irreducible simple subalgebras of su(7), where the 
corresponding irreducible representations are specified by highest wcig hts. 7 Wc 
find the following irreducible simple subalgebras (see Tab. [7]): su(7) with highest 



7 The definition of the highest weight is discussed in Appendix IC.ll 



50 R. Zeier, T. Schulte-Herbriiggcn 

Table 7. Highest Weights of the Irreducible Representations up to Dimension 16 



dim algebra type highest weight(s) 



2 


SU(2) 


s 


(1) 


3 


au(3) 


u 


(1,0), (0,1) 




su(2) 


o 


(2) 


4 


su(4) 


u 


(1,0,0), (0,0,1) 




8(3(2) 


s 


(1,0) 




3u(2) 


s 


(3) 


5 


au(5) 


u 


(1,0,0,0), (0,0,0,1) 




so(5) 


o 


(1,0) 




su(2) 


o 


(4) 


6 


su(6) 


u 


(1,0,0,0,0), 








(0,0,0,0,1) 




sp(3) 


s 


(1,0,0) 




su(2) 


s 


(5) 




so(6) 


o 


(1,0,0) 




au(3) 


u 


(2,0), (0,2) 


7 


au(7) 




Cl 0) 








(0,0,0,0,0,1) 




80 (7) 


o 


(1,0,0) 




92 


o 


(1,0) 




su(2) 


o 


(6) 


8 


su(8) 


u 


(1,0,0,0,0,0,0), 








(0,0,0,0,0,0,1) 




sp(4) 


s 


(1,0,0,0) 




0u(2) 


s 


(7) 




So(8) 


o 


(1,0,0,0), (0,0,1,0) 








(0,0,0,1) 




su(3) 


o 


(1,1) 




so(7) 


o 


(0,0,1) 


9 


Sll(9) 


u 


(1,0,0,0,0,0,0,0), 








(0,0,0,0,0,0,0,1) 




so (9) 


o 


(1,0,0,0) 




su(2) 


o 


(8) 


10 


su(10) 


u 


(1,0,0,0,0,0,0,0,0), 








(0,0,0,0,0,0,0,0,1) 




sp(5) 


s 


(1,0,0,0,0) 




su(2) 


s 


(9) 




so(10) 


o 


(1,0,0,0,0) 




so(5) 


o 


(0,2) 




su(3) 


u 


(3,0), (0,3) 




su(4) 


u 


(2,0,0), (0,0,2) 




su(5) 


u 


(0,1,0,0), (0,0,1,0) 



dim algebra type highest weight (s) 



11 


au(ll) 


u 


(1,0,..., 0), 








(0, ... ,0, 1) 




so(ll) 


o 


(1,0,0,0,0) 




8U(2) 


o 


(10) 


12 


su(12) 


u 


(1,0,..., 0), 








(0,...,0,1) 




sp(6) 


s 


(1,0,0,0,0,0) 




8U(2) 


s 


(11) 




80(12) 


o 


(1,0,0,0,0,0) 


13 


bu(13) 




(1, 0, . . . , 0), 








(0,.'..,0,1) 




so(13) 


o 


(1,0,0,0,0,0) 




8U(2) 


o 


(12) 


11 


bu(14) 


u 


(1,0,..., 0), 








(0, ... ,0, 1) 




SP(7) 


s 


(1,0,0,0,0,0,0) 




su(2) 


s 


(13) 




sp(3) 


s 


(0,0, 1) 




so(14) 


o 


(1,0,0,0,0,0,0) 




so (5) 


o 


(2,0) 




<;n^3~) 




(0, 1, 0) 






o 


(oil) 


15 


su(15) 


u 


(1,0,..., 0), 








(0,...,0,1) 




so(15) 


o 


(1,0,0,0,0,0,0) 




su(2) 


o 


(14) 




so(6) 


o 


(0, 1, 1) 




su(3) 


u 


(4,0), (0,4) 




su(5) 


u 


(2,0,0,0), 








(0,0,0,2) 




su(6) 


u 


(0,1,0,0,0), 








(0,0,0,1,0) 




su(3) 


u 


(2,1), (1,2) 


16 


su(16) 


u 


(1,0,..., 0), 








(o, ...,0, 1) 




sp(8) 


s 


(1,0,0,0,0,0,0,0) 




su(2) 


s 


(15) 




sp(2) 


s 


(1,1) 




so(16) 


o 


(1,0,0,0,0,0,0,0) 




so (9) 


o 


(0,0,0,1) 




so(10) 


u 


(0,0,0,1,0), 








(0,0,0,0, 1) 
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weights 8 (1,0,0,0,0,0) and (0, 0, 0, 0, 0, 1), so (7) with (1,0,0), g 2 with (1,0), as 
well as su(2) with (6). We conclude from Theorem l49l that the given irreducible 
representations of su(7) are unitary and that all the other ones are orthogonal 
(see Tab. [7]). It follows that so(7) is directly embedded in su(7): 



The algebras 02 and su(2) are embedded in so (7), but we still have to determine 
the inclusion relations. All algebras not listed with the corresponding highest 
weight in Tab.|6]are directly contained cither in su(N), sp(N/2) [for ./V even], or 
in so (A) depending on whether the irreducible representation is unitary, sym- 
plectic, or orthogonal. We find the algebra su(2) with the highest weight (6) in 
the third row of Tab. Thus the algebra su(2) is contained in so (7) but only 
indirectly so — via rj 2 : 



Example 51. Consider the case of N — 16. First, we obtain all the irreducible 
simple subalgebras of su(16) (see Tab. [7]): su(16) with highest weights (1,0,..., 0) 
and (0, . . . , 0, 1) as well as so(10) with (0, 0, 0, 1, 0). The cases of irreducible sym- 
plectic representations are sp(8) with highest weight (1,0,0,0,0,0,0,0), su(2) 
with (15), and sp(2) with (1, 1). The irreducible orthogonal representations are 
given by so(16) with highest weight (1, 0, 0, 0, 0, 0, 0, 0) and so(9) with (0, 0, 0, 1). 
We immediately conclude that the algebras sp(8), so(16), and so(10) are directly 
embedded in su(16): 



It follows that all the other cases are directly contained either in sp(8) or so(16) 
as they are not listed in Tab. [S] 



D. Alternating and Symmetric Squares of Representations 

In this Appendix, we enumerate for all compact semisimplc Lie algebras those 
representations whose alternating and symmetric squares are both irreducible. 
We obtain that only the standard representation of su(£ + 1) with I > has 
this property This result is used in Sec. 17.41 We freely use the notation of 
Appendix ICl 

Assume that </> is a representation of a compact semisimplc Lie algebra q on a 
finite-dimensional vector space V with basis {ui, . . . ,Vk}. The representation <f> 
is given as a map from g to the set of complex k x k matrices gl(k, C). Starting 

8 The two irreducible representations of su(7) are related by an inner automorphism. 



■ ■ ■ <-so(7)<-su(7) 



su(2)<-g 2 «-so(7)<-su(7) 



■ ■ ■ <-so(16)<--su(16) 
so(10)* 
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from the representation cf> we construct its tensor square 4>® 2 = <f> ® Ife + ® (f> 
which acts on the fc 2 -dimensional vector space V<£> V with basis {v{ 1 ©v,; 2 | *i , *a G 
{1, . . . , fc}}. This action is defined on the basis by (g € g) 

0® 2 (ff)IVu ® : = [<KsMi] ® ^2 + ® [0(fl)u< a ] 

and it can be extended to the full vector space V © V by linearity. Now we can 
define for <f> its alternating square Alt 2 := </>® 2 |Ait 2 v by restricting </>® 2 to the 
k(k — l)/2-dimensional subspace Alt V C V © V with basis 

{«»! <& v i2 - v i2 © | ii,i 2 G {1, . . . , fc} and ii 7^ i 2 }. 

It is clear that Alt 2 is well defined as (g £ 0) 

(AltVXfiOfaii ® v 42 - w l2 ® UtJ = ® v l2 - v l2 ® [</>(ff)Vu]) 

+ ® [0(s)ui 2 ] - [<£(s)u< a ] ® 

is contained in Alt 2 ^. Similarily, one defines the symmetric square Sym 2 := 
</>® 2 | Sym 2y as the restriction to the fc(fc + I)/2-dimensional subspace Sym 2 V C 
V ®V with basis {v^ © v i2 + v i2 8^ | *i , *2 G {1, . . . , fc}}. We obtain that the 
tensor square <fi® 2 = Alt 2 © Sym 2 decomposes in a direct sum, exactly as the 
tensor product V © V = Alt 2 ^ © Sym 2 F. Dynkin [55] classified the cases when 
Alt 2 is irreducible: 

Theorem 52 (Dynkin). Assume $ is a (finite- dimensional) representation of 
a compact semisimple Lie algebra g. The representation Alt <f> is irreducible if 
and only if <f> is irreducible and the pair (q, (j)) is (up to an outer automorphism 
of q) given in the following table: 



dim(0) Alt 2 dim(Alt 2 0) 



(la) (1,0,..., 0) 2i+l (0,1,0,...,0) (2l+l)i 

(lb) 00(5) (1,0) 5 (0,2) 10 

(2a) (1,0,..., 0) It (0,1,0,..., 0) {2l-\)l 

(2b) so(6) (1,0,0) 6 (0,1,1) 15 

(3) m ^ (0,1,0,...,0) (1,0,1,0,...,0) 3( e + 2 ) 

(4) -(g 1 ) (2,0,...,0) i£±iK£±2) (2,1,0,...,0) 3(*+ 3 ) 

(5) so(10) (0,0,0,1,0) 16 (0,0,1,0,0) 120 

(6) cc (1,0,0,0,0,0) 27 (0,0,1,0,0,0) 351 

(7) (1,0,..., 0) e+i (o,i,o,...,o) 



Proof. If 4> = (f>i © 02 is not irreducible then neither is Alt 2 (</>i © fa) = Alt 2 0i © 
(4>i © 02) © Alt 2 (/> 2 irreducible. The Theorem follows from Thm. 4.7 and Tab. 6 
of Ref. [55]. ■ 



Relying on Theorem 1521 we obtain: 
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Theorem 53. Assume cf> is a (finite- dimensional) representation of a compact 
semisimple Lie algebra g. The representations Alt <f> and Sym 2 </> are both irre- 
ducible if and only if Q = su(£+ 1) with £ > 1 and is (up to an outer automor- 
phism of q) the standard representation [i.e. its highest weight is (1, 0, ... , 0)/. 

Proof. We go through the cases of Theorem [521 Let us denote by 4> x the repre- 
sentation with highest weight x. In the cases (la)-(2b), it follows from Ex. 19.21 

of Ref. |109j that Sym 2 ^^....^) = 0(2, o o) © 0(o,...,o) decomposes. In the case 

of (3), we can use a Pieri-type formula (see Prop. 15.25(h) of Ref. |109j ) to show 
that Sym 2 (O!liO ) = (o ,o,o) © 0(0,2,0) and Sym 2 (O!liU _. iO) = 0( O ,o,o,i,o,...,o) © 
0(o,2,o,...,o) decompose. In the case of (4), we can use again a Pieri-type formula 
(see Prop. 15.25(i) of Ref. [109] ) to show that Sym 2 ^>( 2 ,o,...,o) = 0(o,2,o,...,o) © 
0(4,0,. ..,o) decomposes. In the cases (5) and (6), we explicitly compute the de- 
composition using computer algebra systems such as Li E [5S] and MAGMA [57] , 
We get for (5) that Sym 2 0( OAO ,i,o) = 0(o, 0,0,2,0) © 0(i,o,o,o,o) and for (6) that 
Sym 2 (liO ,o,o,o,o) = 0(o,o,o,o,o,i) © 0(2,o, 0,0,0,0) ■ In t he case of (7), we use again a 
Pieri-type formula (see Prop. 15.25(i) of Ref. [109] ) to show that (0(i iO o))® 2 = 

0(2,o,. ..,0) ©0(0,1,0,. ..,0) = Sym 2 (1A ._ o) ©Alt 2 (1A ._ o) . Therefore, case (7) is the 
only case for which both the alternating and symmetric square are irreducible. 
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